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ABSTRACT 


The  research  described  encompasses  six  tasks  each  of  which  treats,  from 
either  an  experimental  or  theoretical  base,  the  measurement  or  interpretation 
of  the  high- temperature  thermophyslcal  behavior  of  model  systems  which  show 
potential  applications  to  the  environmental  resistant  materials  and  the 
spacecraft  survivability  programs.  While  there  is  significant  overlap 
between  these  tasks,  the  first  two  chapters  are  addressed  to  the 
environmental  resistant  materials  programs,  while  the  later  four  chapters  are 
addressed  to  the  spacecraft  survivability  program. 

Measurements  have  been  initiated,  in  a  specially  designed  apparatus 
using  state  of  the  art  pyrometery,  to  measure  the  triple  point  of  graphite  by 
a  pulse  heating  method.  Preliminary  results  indicate  that  the  surface 
radiance  temperature  of  melting  graphite  at  effective  wavelengths  between  0.5 
and  0.9  micrometers  are  in  the  range  of  4300  to  4400  K. 

Thermodynamic  properties  of  silicon  nitride  (o,  0)  and  boron  nitride 
(hex,  cub)  have  been  determined  to  1300  K.  The  heat  capacity  of  bcc  calcium 
has  been  determined  In  the  range  of  273  to  929  K.  The  enthalpy  of  the  phase 
transition  to  fee  structure  has  also  been  determined.  These  data  are  to  be 
used  to  test  theories  and  models  for  predicting  the  thermodynamic  properties 
of  metals. 

Calculatlonal  algorithms  have  been  developed  to  calculate 
photodiasociation  cross  sections  for  broad  band  continuum  transitions, 
including  molecular  dissociation.  Preliminary  calculations  have  been 
performed  to  study  the  pressure  broadening  of  atomic  and  molecular  species  in 
strong  laser  fields.  Franck-Condon  matrix  elements  for  bound- continuum 
transitions  are  being  used  to  obtain  estimates  of  total  cross-sections. 


i  ii 


A  combined  theoretical  and  experimental  [X'ogram  Is  being  pursued  to 
Increase  the  understanding  of  the  relationship  of  structure  and  reactivity 
of  graphitic  structures.  Kinetics  of  the  phenyl  radical  deposition  on 
pyrolytic  graphite  surfaces  have  been  investigated.  Results  of  this 
investigation  have  been  Interpreted  In  terms  of  the  number  of  available 
reactive  sites. 

The  rates  and  mechanism  of  Ionic  polymerization  In  highly  unsaturated 
organic  compounds,  structures  and  thermochemistry  of  relevant  Ions, 
photodlssociatlon  of  the  polymeric  Ions,  and  the  repair  mechanisms  which 
generate  polymeric  Ions  after  photodlssociatlon  have  been  studied.  The  work 
focusses  on  the  elucidation  of  the  gas  phase  polymerization  and  reactions 
mechanisms  following  laser  Irradiation  Induced  breakdown  In  aromatic 
hydrocarbons,  particularly  those  with  acetylenic  and  oleflnic  side  chains 
such  as  phenylacetylene  and  phenylethylene  (styrene).  Kinetics  and 
mechanisms  of  the  reactions  of  phenyl  Ions  were  also  examined,  and  structures 
of  this  species  generated  from  a  number  of  organic  molecules  have  been 
elucidated. 

Development  of  apparatus  and  techniques  to  study  the  molecular  basis 
for  laser  Induced  vaporization  of  refractory  materials  such  as  graphite 
continued.  Special  features  in  the  vaporization  of  graphite  have  been 
examined.  The  gas  phase  species  from  zo  C5  have  been  observed.  The 
computation  of  thermodynamic  temperature  using  the  ratios  of  these  species 
observed  using  mass  spectrometrlc  techniques  have  been  elucidated.  The 
effects  of  laser  energy  on  the  mass  spectrum  have  also  been  studied. 


A.  Cttzalrllyan*  A.  P.  Nllllar  and  J.  L.  McClura 

AggTPA<?I 

Haaauraaanta  hava  baan  initlatad  to  accurataly  dataraina 
tha  tripla-point  tanparatura  of  graphita  by  a  rapid  pulaa- 
haating  tachniqua.  Tha  baaic  aathod  invoivaa  haating  tha 
apaciaan  in  a  praaaurizad  inart  gaa  anvironaant  (14  MPa>  froa 
rooa  taaparatura  to  tha  aalting  taaparatura  in  laaa  than  0.1  a 
by  tha  paaaaga  of  an  alactrical  currant  pulaa  through  it;  and 
aiaultanaoualy  dataraining  tha  apaciaan  taaparatura  avary  0.1  aa 
by  aaana  of  a  nawly-conatructad  pyroaatar  which  aaaauraa  tha 
apactral  radianca  taaparaturaa  in  aix  wavabanda.  Final  taating 
and  calibration  of  tha  aulti-wavalangth  pyroaatar  ia  naar 
coaplation.  Praliainary  raaulta  indicata  that  tha  aurfaca 
radianca  taaparaturaa  of  aalting  graphita  at  affactiva 
wavalangtha  batwaan  O.S  and  0.9  pa  ara  approxiaataly  in  tha 
ranga  4300  to  4400  K. 

1.  IWTRODUCTIOW 

Tha  aolid-liquid-vapor  tripla  point  of  graphita  has  baan 
tha  aubjact  of  nuaaroua  atudiaa,  particularly  during  tha  paat 
dacada  or  ao;  raviawa  of  tha  invaatigationa  may  ba  found  in  tha 
litaratura  Cl. 21.  Although  aoat  aaaauraaanta  yiald  a  tripla- 
point  praaaura  of  about  110  ata.  conaidarabla  controvaray  atill 
axlata  aa  to  tha  tripla-point  taaparatura  with  raportad  vaiuaa 
apanning  tha  ranga  of  approxiaataly  4000  to  5000  X.  A  aa^or 
problaa  in  auch  aaaauraaanta  ariaaa  froa  tha  high  vapor  praaaura 
of  graphita:  a  pluaa  of  coolar  carbon  vapor  (in  tha  fora  of 
aoot>  foraa  naar  tha  apaciaan  aurfaca  aa  tha  tripla  point  la 
approachad.  tharaby  aavaraly  liaiting  tha  accuracy  of  optical 
pyroaatry . 

A  faaaibility  atudy  C33  waa  parforaad  racantly  in  our 


laboratory  to  dataraina  if  tha  af facta  of  apaciaan  avaporation 


on  pyronotric  taaporatur*  moaauroaont.a  could  ba  auppraaaad,  or 


at  laaat  nlnlaizad.  Surface  radiance  data,  obtained  by  aeana  of 
a  alngle-wavelength  (0.6S  pa)  hlgh-apeed  pyroaeter,  ahowed  that 
a  reaaonably  well-defined  aelting  plateau  in  the  teaperature- 
tine  function  could  be  obtained  under  certain  conditiona:  (1)  a 
anall  anount  of  oxygen  <10  to  30a)  aust  be  added  to  the  inert 
gaa  environaent  of  the  apecinen*  The  carbon  vapor  conbinaa  with 
the  oxygen  to  fora  optically  tranaparent  gaaea  CO  and  CO2; 

(2)  the  apeciaena  nuat  be  rapidly  heated  at  ratea  of  about 
103  K^a^^or  higher.  The  very  ahort  heating  tiaea 
(leaa  than  0.1  a)  tend  to  nininize  the  effecta  of  convective 
turbulence  in  the  aighting  path  of  the  pyroaeter. 

Aa  a  raault  of  thia  atudy,  a  aulti -wavelength  pyroaeter 
waa  deaigned  and  conatructed  C4]  to  provide  a  aeana  of 
eatiaating  the  true  tenperature  of  a  apeciaan  froa  aurface 
radiance  data.  The  final  teating  and  calibration  of  thia 
inatrunent  ia  nearing  coapletion. 

In  the  preaent  report,  we  briefly  deacribe  <1>  the  aa3or 
conponenta  of  the  aeaaurenent  ayatea  and  (2)  preliainary  reaulta 
obtained  froa  aeaaureaenta  on  graphite  with  the  aulti -wavelength 
pyroaeter  at  teaperaturea  in  the  vicinity  of  the  triple-point. 

2.  MEASUREHEWT  SYSTEM 

The  high-apeed  aaaaureaent  ayatea  includea  an  electric 
power-pulaing  circuit  and  ita  aaaociated  aeaauring  and  control 
eireuita,  a  high-praaaura  experiaent  chaaber,  a  aulti-wavelangth 
pyroaeter,  and  a  data  acquiaition  ayatea.  A  functional  diagram 


of  tho  ■•aauromont.  ayataii  !•  ahown  in  Fig.  1 
2.1.  Pulainq  Circuit 


Tha  powar-pulaing  circuit  conaiata  of  tha  apaciaan  in 
aariaa  with  a  battary  bank,  an  adjuatabla  raaiatanca  (watar- 
coolad  Inconal  tuba),  and  a  faat  acting  awitch.  Adjuataanta  to 
tha  battary  bank  voltaga  <up  to  38  V>  and  to  tha  langth  <hanca. 
raaiatanca)  of  tha  Inconal  tuba  in  tha  circuit,  prior  to  pulaa 
haating.  anabla  control  of  tha  apaciaan  haating  rata.  Tha 
tiaing  of  varloua  avanta.  auch  aa  cloaing  and  opaning  of  tha 
awitch  and  triggering  of  alactronic  Inatruaanta.  ia  achiavad  by 
aaana  of  tiaa-dalay  pulaa  ganaratora. 

2.2.  Exoarlaant  Chaabar 

Tha  axpariaant  chaabar  la  an  ultra-high  praaaura  call 
capable  of  providing  an  inert  gaa  anvironaant  at  praaauraa  up  to 
200  HPa  ( 2000  ata) .  Two  optical  porta  through  tha  cylindrical 
wall  on  oppoaita  aidaa  of  tha  vaaaal  (viewing  cona  of  about 
200  )  facilitate  tha  alignaant  of  tha  apaciaana  and  provide 
a  aighting  path  for  tha  pyroaatar  through  aapphira  windowa 
(diaaatar.  25  aa;  thicknaaa.  32  aa) .  Tha  apaciaan/aount 
aaaaably  ia  aupportad  froa  tha  cloaura  head  which  containa 
alactrlcal  faadthrougha  for  tha  pulaa-haating  currant. 

Tha  apaciaan  ia  claapad  onto  two  alactrodaa  in  tha 
aounting  aaaaably:  a  atationary  lower  electrode  and  an  upper 
electrode  which  ia  connactad  to  an  axpanaion  joint  (flaxibla 
phoaphor-bronza  ballowa).  Thia  arrangaaant  allowa  for  tharaal 
axpanaion  of  tha  apaciaan  along  ita  axial  direction  during  rapid 
pulaa  haating. 


2.3.  SD«ci««n« 

Th«  •p«cia«n  aat«ri«l,  dMlgnatad  by  tha  aanufacturar  aa 
POCO  AXM-5Q1  graphita*  waa  auppllad  in  tha  fora  of  thin  <0.5  aa) 
ahaata.  Tha  apaciaana  wara  fabricatad  into  raetangular  atripa 
with  noainal  diaanaiona  of  25  aa  long  by  3  aa  wida.  Tha  cantar 
portion  of  aach  apaciaan  waa  than  **nackad-down**  to  dafina  an 
“affactiva**  portion  of  tha  apaciaan  which  undargoaa  aalting. 

Tha  diaanaiona  of  tha  affactiva  apaciaan  wara  noainal ly:  langth, 
2  aa;  width,  1  aa;  thlchnaaa,  0.5  aa. 

2.4.  Multi -Wavalancth  Pvroaatar 

Tha  pyroaatar  ia  capabla  of  aaaauring  tha  apactral 
radianca  taaparaturaa  of  a  rapidly  haating  apaseiaan  in  aix 
wavabanda  with  a  tlaa  raaponaa  of  50  ;jia  (for  aach  apactral 
band).  Tharaal  radiation  froa  a  circular  targat  araa  (0.5  aa 
diaaatar)  on  tha  apaciaan  ia  focuaaad  by  tha  lana  ayataa  in  tha 
pyroaatar  onto  a  fibar  optic  bundla.  Tha  lattar  ia  dividad 
(noainally  at  randoa)  into  aix  aaparata  bundlaa  which  illuainata 
aix  intarfaranca  filtara  and  ailicon  dioda  datactora.  Fiva  of 
tha  intarfaranca  filtara  hava  a  apactral  bandwidth  of  about 
0.1  pa,  cantarad  at  0.5,  0.6,  0.7,  0.8  and  0.9  pa.  Tha  aixth 
filtar  haa  a  bandwidth  of  0.03  pa  cantarad  at  0.65  pa.  Tha 
photocurranta  froa  tha  diodaa  ara  convartad  by  high-atability 
aaplifiara  to  proportionata  voltagaa,  which  ara  awitchad 
aaquantially  avary  0.1  aa  through  a  aultiplaxar  to  tha  pyroaatar 
output. 


Tha  calibration  of  tha  pyroaatar  involvaa  two  atapa. 
Firat.  tha  0.65  pa  channal  ia  calibratad  in  a  ataady  atata 


•xp*riM«nt  with  a  tungstan-fllaaant  standard  lamp  which*  in 
turn*  had  baan  calibratad  against  tha  MBS  Photoalactric 
Pyrosatar  oparating  at  tha  sasa  a£factiva  wavalangth.  Tha 
calibration  is  than  transfarrad  to  tha  othar  fiva  channals  by 
dynamically  racording  tha  radiation  from  a  blackbody  hola  in  a 
graphita  tuba  as  it  is  haatad  by  a  pulsa  of  alactric  currant. 

Additional  data! la  concarning  tha  dasign*  oparation  and 
calibration  of  tha  sulti -wavalangth  pyrosatar  ara  givan 
alsawhara  C4] . 

2.5.  Data  Acquisition 

Tha  nultiplaxad  analog  signals  from  tha  pyrosatar  ara 
digitizad  (IS-bit  rasolution)  and  storad  by  a  digital 
oscilloscopa  having  a  sasory  capacity  of  16X  data  points.  Tha 
tisabasa  of  tha  digital  oscilloscopa  is  suppliad  by  tha 
pyroaatar  alactronics  which  provida  a  saquanca  of  TTL-cospatibla 
"sasory-addrass  advanca"  pulsas  synchronous  with  tha  sultiplaxer 
switching.  After  each  axparisant*  tha  data  is  transfarrad  from 
the  oscilloscopa  to  an  off-line  computer  for  subsequent 
analyses . 

3.  MEASUREMENTS  AND  PRELIMINARY  RESULTS 

During  tha  course  of  successive  calibration  exparimants. 
tha  new  pyrosatar  was  found  to  exhibit  a  long-tars  drift  in  its 
response.  As  a  result  of  extensive  tasting  of  pyrosatar 
cosponants*  tha  source  of  instability  was  detersinad  to  ba  tha 
intarfaranca  filters.  Raplacasant  filters  ware  ordered  and 


should  ba  available  early  in  FY86.  In  spita  of  tha  instability 


how«v«r,  um£u1  inforaation  could  atlll  ba  obtainad  with  tha 
pyroaatar  and  tharafora,  tha  aarlaa  of  pulaa  axpariaanta  to 
dataraina  tha  trlpla-point  taaparatura  of  graphita  waa  bagun. 

Prior  to  aach  pulaa  axpariaant*  a  atrip  apaciaan  waa 
aountad  in  tha  axpariaant  call  which  waa  than  praaaurizad  to 
14  MPa  ( ~  140  ata)  with  an  argon/oxygan  gaa  aixtura;  tha  oxygan 
contant  variad  with  axpariaant  batwaan  5  and  33  x.  Tha  apaciaan 
waa  than  rapidly  haatad  froa  rooa  taaparatura  to  ita  aalting 
taaparatura  in  laaa  than  0.1  a  by  paaaing  an  alactrical  currant 
pulaa  through  It.  Tha  duration  of  tha  currant  pulaa  rangad  froa 
approxiaataly  40  to  70  aa.  Haating  rataa  wara  noainally  batwaan 
l.S  and  3  x  105  K*a-1. 

Tha  raaulta  of  a  typical  pulaa  axpariaant  ara  praaantad  in 
Fig.  2  which  illuatrataa  tha  variation  of  apactral  radianca  froa 
tha  apaciaan  aurfaca  with  tiaa,  aa  aaaaurad  by  tha  aix  channala 
of  tha  pyroaatar.  Plataaua  in  tha  radianca  varaua  tiaa 
functiona  Indlcata  that  tha  apaciaan  haa  undargona  aalting. 
Further  evidanca  of  aalting  aay  ba  aaan  in  Fig.  3  which  ahowa 
that  tha  aurfaca  of  tha  raaaining  portion  of  affactiva  apaciaan 
la  glaaay  in  appaaranca. 

For  a  given  axpariaant,  radianca  taaparaturaa  along  tha 
flat  portion  of  tha  aalting  plataau  wara  avaragad  for  aach 
channal .  A  typical  aaan  value  Involved  approxiaataly  15 
taaparaturaa  and  a  standard  deviation  of  about  10  K. 

Tha  naan  values  of  radianca  taaparatura  corasponding  to 
tha  0.&5  >ia  channal  ara  plotted  In  Fig.  4  as  a  function  of 
oxygan  contant  In  tha  prasaurizlng  gas  anvlronnant.  As  nay  ba 


•Mil,  thm  mMsurad  radionca  taaparafcura  (at  0.6S  >im>  o£  graphita 
at  tha  trlpla  point  daeraaaaa  linaarly  with  daeraaaing  oxygan 
contant.  A  ainilar  bahavior  waa  obaarvad  in  tha  raaulta 
racordad  by  tha  othar  fiva  ehannala.  Thia  affact  appaara  to 
ariaa  from  tha  raaction  o£  axcaaa  oxygan  with  carbon  in  tha 
apaciaan  aurfaca,  which  auggaata  that  tha  "corract**  valua  of 
radianca  taaparatura  at  aalting  aay  ba  datarainad  by 
axtrapoiation  of  tha  aaaaurad  raaulta  to  zaro  oxygan  contant. 

Tha  praliainary  raaulta  indicata  that  tha  radianca 
taaparaturaa  of  graphita  at  ita  tripla  point  for  affactiva 
wavalangtha  batwaan  O.S  and  0.9  >ia  ara  approxiaataly  in  tha 
ranga  4300  to  4400  K«  Howavar,  no  raliabla  aatiaata  of  tha  trua 
tripla-point  taaparatura  could  ba  aada  bacauaa  of  tha 
uncartaintiaa  craatad  by  tha  unatabla  intarfaranca  filtara  in 
tha  pyroaatar. 

Tha  final  tasting  and  calibration  of  tha  pyroaatar  will  ba 
coaplatad  aarly  in  FYSG  aftar  raplacing  tha  filtara  with  aora 
atabla  units.  Subsaquant  pulaa  axpariaanta  should  than  anabla 
an  accurata  dataraination  of  tha  tripla-point  taaparatura  of 


graphita 


(13  F.  P.  Bundy,  J.  Gnophynicnl  Rnsnarch  BS.  6930  <1960>. 
C2]  H.  A.  Shalndlin,  High  Tanparatura  19.  467  (1981>. 

(33  A.  Cazalrllyan,  A.  P.  Mlillar  and  J.  L.  MeClura,  Raport. 
AFOSZR-ISSA-83-00038,  1983.  p.l. 

(43  G.  H.  Folay,  to  ba  publlahad. 
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Fig.  Oscllloaeopg  tr«e«  photograph  ahoMing  tha  variation 
with  tlaa  of  tha  apactral  radianea  from  tha  apaciaan  aurfaca*  aa 
aaan  by  alx  channala  of  tha  aultl*wavalangth  pyroaatar»  during  a 
typical  axpariaant  in  which  tha  apaciaan  waa  pulaa  haatad  to* ita 
aalting  point  in  a  praaaurizad  gaa  anvironaant  (80h  argon ^  20h 
oxygan)  at  14  HPa  and  at  a  haating  rata  of  3xlO^K>a’‘l. 


1  3 1  A  photograph  of  tha  lowar  portion  of  tha  apaciaan 

waa  pulaa  haatad  in  a  SOX  argon/  20X  oxygan  anvironaant. 
Tha  glaaay  appaaranca  axhibitad  by  tha  aurfaea  of  tha  raaaining 
portion  of  affactiva  apaciaan  indicataa  that  aalting  haa  occurrad 
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HI  CM  TEMPERATURE  ENTHALPY  MEASUREMENTS 
David  A.  Dltmars,  Principal  Investigator 

ty 

I.  Thermodynamic  Properties  of  Environmentally  Resistant 
Ceramics  and  Refractories 


Introduction 

Heat-capacity  calorimetry  is  one  of  most  cost-effective  approaches 
(often,  it  is  the  only  available  approach)  for  obtaining  accurate  fundamental 
thermodynamic  property  data  for  materials  which  are  candidates  for 
applications  involving  high  temperatures  and  corrosive,  hostile  operating 
environments.  Such  data  may  have  direct  engineering  applications  in 
optimizing  heat  transfer,  guarding  against  thermal  shock,  and  in  evaluating 
the  accuracy  of  some  thermophysical  property  data  measured  by  independent 
transient  methods.  It  is  also  used  in  the  calculation  of  high- temperature 
chemical  equilibria  and  in  the  generation  and  analysis  of  ceramic  and  alloy 
phase  diagrams. 

The  NBS  high-temperature  calorimetry  facility,  maintains  an  isothermal 
drop  calorimeter  for  the  temperature  range  273-1200  K  and  an  adiabatic 
receiving  calorimeter  for  the  temperature  range  1200-2600  K.  The  temperature 
range  of  the  latter  instrument  overlaps  that  of  the  NBS  millisecond  pulse 
calorimeter  and  provides  a  valuable  check  on  the  transient  measurement 
approach  of  the  pulse  technique.  Substances  investigated  may  be  liquid  or 
solid  and  include  both  electrical  conductors  and  non-conductors.  MBS  has 
employed  this  facility  to  investigate  refractory  and  other  materials  of 
interest  to  the  Air  Force.  Attempts  are  made  to  identify  materials  for  which 
currently  available  thermodynamic  data  are  Incomplete  or  unreliable  due  to 
poor  characterization  of  the  samples  or  inaccuracies  in  measurement  • 
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technique. 
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Nitrides,  borides,  and  carbides  form  three  very  Important  classes  of 
refractory  materials  finding  application  in  high-temperature  design  and  as 
components  of  composite  materials.  Silicon  nitride  and  boron  nitride,  for 
instance,  are  two  Important  refractory  construction  materials  for  which  the 
thermodynamic  data  is  of  questionable  reliability.  Heat  measurements  on 
these  materials  are  sparse  and  their  Interpretation  is  complicated  by  '^he 
simultaneous  occurrence  of  two  crystalline  forms  of  each  of  these  substances. 

Measurements 

The  enthalpy  of  silicon  nitride  (a, 8)  relative  to  273  K  and  that  of 
boron  nitride  (hex,  cubic)  relative  to  273  K  have  both  been  measured  to 
1200  K.  The  two  crystalline  forms  of  Si3M4  are  based  on  a  hexagonal  unit 
cell  and  differ  principally  in  having  a  slightly  different  arrangement  of  the 
nitrogen  tetrahedra  enclosing  each  Si  atom.  The  8-phase  is  thought  to  be 
metastable  at  room  temperature.  Mo  low  temperature  (below  273  K)  heat 
capacity  data  on  812^4  have  appeared  in  the  literature.  Only  two  high 
temperature  heat  capacity  studies  for  SI2N4  have  been  reported  [i,2]  and 
these  differ  by  as  much  as  30)1  below  13OO  K.  No  structural  characterization 
of  the  material  Investigated  in  [1,2]  was  given.  The  present  samples  of 
aiSi2N4  and  8SI2N4  were  submitted  for  both  structural  and  chemical  analysis. 
They  showed  98“^  percent  chemical  purity  and  i  95  percent  phase  purity.  The 
heat-capacity  data  derived  from  the  experiments  are  presented  below: 

aSi^N4  = 

Cp  -  (-0.3262989)x10“'*t2  +  (0.6552924)x10“’t  ♦  (0. 1  38441 1  )x1 0^ 

*  (-0.1558073)x105  x  (273.15/t2)  (1) 

3S1:^N4: 

Cp  -  (-0.5l83574)x10“^t2  +  ( 0 . 9274656)  1 0"H  ♦  (0. 1 295339)x1  o3 
+  (-0.1324359)x105  x  (273.15/t2) 


(2) 


where 


Cp  -  heat  capaclty/J*K“^ 
t  -  Celsius  temperature 
T  -  Kelvin  temperature  -  t  ♦  273.15 

The  two  crystalline  forma  of  boron  nitride  are  a  hexagonal  ("graphitic") 
phase  and  a  cubic  ("tetragonal”,  or  "dlamond-llke" )  phase.  The  stable  phase 
at  atmospheric  pressure  and  temperatures  above  500  K  Is  the  hexagonal  phase. 
The  situation  with  BN  Is  analogous  to  that  which  exists  with  carbon 
(graphlte/dlamond) .  Cubic  BN  Is  produced  by  shock  compression  at  high 
temperatures  and  pressures  and,  like  diamond.  Is  stable  at  a  temperature  of  a 
few  hundred  Kelvin  and  below.  The  hlgh-temperature  limit  of  stability  of 
cubic  BN  is  uncertain  but  it  appears  likely  to  revert  to  the  hex  form  below 
1700  K.  Fortunately,  there  exist  single  low-temperature  (4-300  K)  studies 
for  both  forma  of  BN  [3# 4]. 

The  chemical  purity  of  both  BN  samples  was  better  than  99"^  (mass) 
percent.  The  heat  capacity  data  derived  from  the  experiment  are  given  below. 
Hexagonal  BN; 

Cp  -  (-0.2385988)x10~'*t2  +  (0.4458208)x10"’t  ^  (0.2511 51 5)x1o2 

>  (-0.2267366)x10‘*x(273.15/t2)  (3) 

Cubic  BN; 

Cp  -  (-0.2509428)x10“'’t2  ♦  (0. 455799«)x1 0" 't  ♦  (0.2637653)x1o2 

+  (-0.4000692)x10^x(273.15/t2)  (4) 

where, 

Cp  -  heat  capacity/J*mol"^ 

t  -  Celsius  temperature 

T  -  Kelvin  temperature  -  t  +  273.15 

The  high-and  low-temperature  (below  room  temperature)  heat  capacity  data  are 


in  good  agreement . 
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II.  Heat  Capacity  of  bcc  Calcium;  Enthalpy  of  fcc/bcc  Transition 


Introduction 

The  existence  of  an  unusual  rise  (up  to  50f  or  more)  in  the 
high-temperature  heat  capacity  of  several  bcc  refractory  metals  has  resisted 
theoretical  Interpretation  to  data  [1].  Attempts  have  been  made  by 
theoreticians  to  define  more  adequate  models  (choice  of  interatomic 
potential,  radius  of  neighbor-neighbor  interactions,  etc.)  and  these  have 
been  tested  on  simple  fee  metals  for  which  there  exists  substantial 
experimental  thermodynamic  data  base  [2].  In  the  case  of  aluminum,  the 
comparison  [3]  of  theoretical  results  for  Cy(3)  to  the  values  for  Cv(s), 
derived  from  the  literature,  cast  doubt  upon  the  accuracy  of  the  literature 
data.  Mew,  definitive  measurements  of  Cp(s)  in  the  range  273  K  -  929  K  [4] 
indicated  the  literature  data  to  be  in  error  and  yielded  a  highly 
satisfactory  confirmation  of  the  theoretical  model  [5].  In  order  to  extend 
the  fee  model  to  bcc  metals,  it  is  desirable  to  obtain  accurate  heat-capacity 
data  on  a  metal  exhibiting  both  fee  and  bcc  phases  at  high  temperatures, 
without  complicating  magnetic  transitions.  Such  behavior  is  found  in  the 
alkaline  earths  calcium  and  strontium.  The  thermodynamic  functions  for  these 
highly-reactlve  metals  are  sparse  and  Inaccurate  due  to  experimental 
difficulties  in  achieving  and  maintaining  adequately  low  sample- impurity 
levels..  We  report  new  data  on  a  pure  calcium  sample. 

Samples 

Triply-distilled  calcium  metal  was  obtained  from  Rare  Earths  Products, 
Ltd.,  Widnes,  Cheshire,  U.K.  This  sample  ("REP")  was  analyzed  by  several 
techniques  for  the  major  impurities.  It  was  then  re-processed  by  a  modified 


distillation  method  at  the  Materials  Preparation  Center  of  the  State 


University  of  Iowa.  The  re-processed  sample  is  referred  to  as  "REP/AMES". 
Table  1  presents  the  measured  impurity  levels  or  the  two  samples.  A  specimen 
of  sample  REP/AMES  of  mass  2.19665  g  was  encapsulated  in  a  tantalum  container 
in  an  argon  atmosphere.  The  tantlalum  container  was  sealed  in  a  Pt/IORh 
container  in  an  inert  atmosphere. 


Table  1 

Calcium  Impurity  Levels 

REP  REP /AMES 


At  X 

Mole  % 

At  X 

Mole  X 

0 

2.606 

2.238 

H 

0.737 

0.1*77 

C 

0.0U7 

m 

N 

0.021 

m 

Ba 

0.012 

m 

CaO 

2.676 

0.239 

CaH2 

2.676 

0.2*10 

Ca 

96.953 

99.521 

Methods 

The  enthalpy,  relative  to  273.15  K,  of  the  sample  plus  encapsulation  was 
measured  in  a  precise,  isothermal,  phase-change  calorimeter  over  the 
temperature  range  673  -  1112  K  using  previously  documented  methods  [5,7]. 
Multiple  measurements  were  made  at  several  temperatures  distributed  equally 
throughout  the  temperature  range  as  a  ahecK  on  the  precision  of  measurement. 
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In  addition,  an  earlier  series  of  enthalpy  measurements  on  an  empty  Pt/lORh 


capsule  nearly  identical  to  the  one  used  to  encapsulate  the  calcium  specimen 
was  available  [8].  The  specimen  enthalpy  data  were  obtained  by  subtracting 
at  each  temperature,  the  measured  specimen  plus  encapsulation  enthalpy,  from 
that  of  the  encapsulation  alone.  The  contribution  of  the  tantalum  to  the 
gross  measured  heat  was  calculated  using  the  tantalum  enthalpy  data  of  the 
JANAF  tables  [9].  Corrections  were  made  to  the  gross  measured  sample  heats 
for  the  small  mass  difference  between  the  specimen  and  empty  capsules  and  for 
the  heat  due  to  the  major  assumed  impurities,  CaO  and  CaH2  (cf  Table  1).  In 
addition,  a  second  calcium  specimen  was  independently  encapsulated  in  a 
tantalum  vessel  and  the  solid-solid  transition  temperature  measured  by 
thermal  analysis. 

The. enthalpy  data  for  the  fee  phase  were  fitted  with  a  linear  model: 

(Ht  -  H273.15)  -  AT  ♦  3  J‘mol"’  (1) 


where , 

H  »  enthalpy 
T  -  Kelvin  temperature 

The  coefficients  and  the  standard  deviation  (oj.)  of  this  fit  are: 

A  -  *23.0788 

3  -  -77H2.9'* 

iJf  =  35.52 

The  enthalpy  data  for  the  bcc  phase  were  fitted  with  a  quadratic  model: 
(H.J.  -  H273,i5)  -  CT^  *  DT  ♦  E  J*mol~’  (2) 

The  coefficients  and  the  standard  deviation  (0^)  of  this  fit  are: 

C  -  *0.376302S-02 
D  -  *0.154362E+02 


Results 


An  fcc-bcc  transition  temperature  of  71 5  K  was  measured.  The  corrected 
net  enthalpy  data  for  the  specimen  are  given  in  Table  2. 


Table  2.  Enthalpy  Data  for  CaCs)  ^ 
K _ (H.^  —  /  J  POl 


683.244 


11484.99 


683.038 


11472.51 


673.019 


11129.12 


672.961 


11125.08 


703.030 


11997.28 


702.900 


11997.28 


702.900 


11978.14 


773.292 

773.244 

822.55 


15003.70 

15018.01 


15485.00 


873.03 


18019.52 


872.89 


17992.36 


923.71 


19590.94 


974.37 


21238.16 


974.28 


21 1 74.28 


1023.86 


22848.29 


1073.89 


24498.26 


1104.44 


25630.67 


M.W.  =  40.08 
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In  Table  3t  the  present  results  have  been  compared  to  selected  values 
obtained  from  the  most  widely  used  compilations  of  thermodynamic  properties 
[9,10,11 ]. 


Table  3-  Summary  of  Selected 
Thermodynamic  Data  for  Calcium 


Ref.  [10] 

Ref.  [11] 

Ref.  [9] 

Present  Work 

^fcc/bcc 

720+  3 

716 

721  ±  2 

715+10  /K 

^^fcc/bcc 

220±50 

222 

220±50 

238±  4  /J-moT 

^^Pfcc/bcc 

-9.4 

-3.8 

+5.3 

-4.3  /  % 

Cp(Tm) 

45.52 

43.20 

41 .21 

34.98  /J*moT 

-W-1 


.  -v 

•  *  *  A  •  • 
.  •  A  • 
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•  -^4 
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Broad-Band  Spectroscopy  o-f  Small  Molecules  in 
Intense  Laser  Fields 


Principal  Investigators:  P.  S.  Julienne  and  F.  H.  Mies 


Introduction 

The  general  goal  of  this  research  is  to  provide  calculations 
of  radiative  cross-sections,  fragment  state  distributions,  and 
intense  laser  saturation  effects  for  broad-band,  wavelength 
insensitive  continuum  spectral  features  of  small  gaseous  molecules 
of  possible  interest  in  spacecraft  survivability. 

Theoretical  input  into  this  program  is  provided  by  extensive 
sets  of  computer  codes  available  to  the  quantum  chemistry  group. 
State  of  the  art  calculations  of  necessary  molecular  potentials, 
couplings  and  radiative  transition  dipoles  are  assessible. 
Scattering  codes  provide  cross-sections  and  branching  ratios  of 
photodi ssoci ati on  fragments.  Programs  for  evaluation  of  adiabatic 
Franck-Condon  matrix  elements  and  non-ad i abat i c  radiati/e 
scattering  matrix  elements  allow  for  evaluation  of  attenuation 
coefficients  for  laser  radiation  propagating  through  gasecus  vapors 

In  addition  to  numerical  procedures,  the  program  requires 
extensive  analysis  of  many  new  theoretical  features  introduced  pv 
verv  intense  laser  fields. 

Progress  Puri na  FY  35 

•;  1  >  A  new  CDC  Cyber  205  central  computer  facility  has  been 
installed  at  NBS  and  a  coed  deal  of  our  effect  has  been  devoted  to 
recompiling  our  cooes  on  this  vastly  improved,  sixth  generation 
computer.  In  addition  to  transferring  e;;isting  programs,  we  are 
exploiting  the  unique  capabilities  of  the  205,  especially  with 
regard  to  vectorising  the  scatte'^ing  codes.  This  upgrading  of 
our  programs  has  profoundly  expanded  our  capabilites  to  treat  the 
Dhctodi  ssoci  ati  on  o-^  ool/atomic  svstems  with  manv  coupled  channels. 
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(2)  Codes  have  been  modified,  and  preliminary  calculations  have  been 
performed  to  study  the  pressure  broadening  of  atomic  and  molecular  lineshapes 
in  extremely  strong  laser  fields.  Although  this  represents  the  spectroscopic 
equivalent  of  a  continuuitt-continuum  transition,  the  half-collision  aspects  of 
the  calculations  are  completely  transferable  to  forthcoming  calculations  of 
bound-continuum  photodlssoclatlon  cross-sections.  Manifestlons  of  a.c.  stark 
shifts  and  the  laser-induced  Mollow  llneshape  associated  with  the  dressing  of 
fragments  by  the  radiation  field  are  extracted  from  the  numerical  results. 
This  suggests  that  normal  perturbation  theory  as  applied  to  radiative 
interactions  must  be  carefully  scrutinized  when  applied  to  laser  attenuation 
mechanisms. 

(3)  The  testing  of  various  algorithms  for  calculation  of  photodissociation 
cross-sections  introduced  into  tne  multichannel  close-coupling  codes  is 
nearly  complete.  We  are  settling  on  a  numerical  procedure  devised  by  Singer 
etal^  which  allows  for  proper  description  of  both  the  total  cross-sections 
and  the  branching  into  fragment  states.  This  gives  very  stable  procedures 
for  evaluating  bound-continuum  transitions,  including  all  possible  subtleties 
associated  with  non-adlabatic  effects  in  molecular  dissociation.  Most  of  the 
important  features  of  non-linear  response  to  intense  laser  fields  can  be 
simulated  by  these  codes,  although  some  well  understood  modifications  might 
have  to  be  introduced  in  future  studies  depending  on  the  radiation  parameters 
employed  in  spacecraft  studies. 

(4)  We  have  applied  our  analytic  treatment  of  molecular  continuum  states  to 
the  analysis  of  molecular  photodlssoclatlon.  This  is  based  on  a  generalized 
version  of  quantum  defect  theory  as  applied  to  predissociation  of 
energetically  excited  molecules.  Given  close  coupled  scattering 
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wavefunctions  obtained  form  numerical  solutions  we  can  extract  analytic 
parameters  which  can  greatly  enhance  our  ability  to  calculate 
photodissociation  cross-sections  over  wide  ranges  of  laser  frequencies. 

(5)  During  the  course  of  our  study  of  strong  field  effects  on 
pressure-broadened  lineshapes,  we  have  developed  a  convergent  analysis  of 
radiative  matrix  elements  which  has  many  applications  throughout  lineshape 
theory. 

(6)  To  complement  our  exact,  fully  non-adiabatlc  evaluation  of 
photodlssoclation  cross-seotiona  for  diatomic  molecules,  we  are  currently 
using  Franck-Condon  matrix  elements  for  bound- continuum  transitions,  and 
theory  previously  developed  for  excimer  transitions^* 3  to  obtain  estimates  of 
total  cross-sections.  Although  these  simple  calculations  should  give  us 
reasonable  estimates  of  laser  attenuation  coefficients  for  diatomic 


molecules,  quantitative  application  to  larger  molecules  remains  to  be  tested. 
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Two  papers  supported  by  this  research  proposal  have  been  submitted  for 
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A  generalized  multichannel  quantum  defect  analysis  of  the  close  coupled 
molecular  wavefunctlon  for  arbitrary  potentials  gives  a  unified  description 
of  the  bound  and  scattering  states  of  the  diatomic  molecule  AB  by  defining  a 
real  symmetric  matrix  T(E)  which  can  be  analytically  continued  across  a 
dissociation  threshold:  The  theory  also  leads  to  a  half  collision  analysis 
of  weak  predissociation  and  photodissociation  whereby  the  transition  ampli¬ 
tude  matrix  la  separated  into  factors  representing  short  range  Franck-Condon 
transitions  and  long  range  final  state  Interactions.  For  sufficiently  high 
fragment  separation  velocity,  the  effects  of  final  state  Interactions  on 
product  distributions  are  well  approximated  by  small  departures  from  the 
recoil  limit  as  prescribed  by  pure  frame  transformations  among  Hund's  cou¬ 
pled  electronic-rotational  states. 


1.  INTRODUCTION 

The  process  with  which  we  are  concerned  here  is  that  of  molecular  dissocia¬ 
tion,  that  is,  the  formation  of  an  unstable  state  of  the  molecule  AB  which 
separates  to  fragments  A  and  B.  There  are  a  variety  of  ways  by  which  a  dis¬ 
sociative  continuum  state  can  be  excited;  through  initial  preparation  of  a 
bound  metastable  state  which  predissoclates ,  through  direct  photoexcitation,  or 
perhaps  through  some  colllslonal  excitation  process.  Generally  the  items  of 
interest  for  measurement  or  calculation  are  the  total  rate  at  which  the  process 
occurs  and  the  branching  ratios  which  describe  the  distribution  of  particular 
final  states. 

The  dissociative  states  of  the  species  AB  are  also  the  collisional  or  scat¬ 
tering  states  of  the  fragments  A  B.  Therefore  there  has  been  a  considerable 
amount  of  interest  and  success  recently  in  generalizing  close  coupling  scatter¬ 
ing  algorithms  for  the  numerical  calculation  of  photodissociation  transition 
amplitudes  (1-7).  The  theory  (8-13)  which  we  will  summarize  here  is  complemen¬ 
tary  to  these  numerical  approaches.  It  is  based  on  an  analysis  of  the  struc¬ 
ture  of  the  close  coupled,  or  multichannel,  time  independent  scattering  wave- 
function.  It  will  hopefully  lead  to  useful  insights  concerning  molecular 
dissociation  aurid  also  lend  itself  to  the  development  of  quantum  and  semi- 
classical  approximations.  Although  our  development  has  concentrated  on  appli¬ 
cation  to  neutral  atom  scattering,  that  is,  diatomic  dissociation,  many  of  the 
features  of  the  theory  are  generail  and  need  not  be  restricted  to  this  case. 

We  call  our  approach  a  generalized  multichannel  quantum  defect  analysis 
(MCQDA)  (10),  since  it  is  strongly  motivated  by  the  concepts  of  standard  multi- 
cnannel  quantum  defect  theory  (MQDT)  (14,15).  The  latter  is  designed  to  treat 
problems  associated  with  bound  Rydberg  states,  autoionizing  resonances,  and 
electron-ion  scattering.  It  has  not  only  been  successfully  applied  to  many 
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atomic  problems,  but  has  also  been  extended  to  treat  molecular  bound  and  dis-  ~~ 
soclative  states  in  the  vicinity  of  an  ionization  threshold  (16,17).  “This 
standard  MQDT  makes  heavy  use  of  the  properties  of  the  Coulomb  potential  and 
its  solutions.  The  generalized  MCQDA  has  been  developed  by  us  (10,11)  and  ''t 

others  (17-19)  to  apply  to  systems  described  by  arbitrary  “potentials;"  for  exam¬ 
ple,  the  Born-Oppenheimer  potentials  for  the  interaction  of  A  and  B. 

There  are  two  features  of  the  theory  which  we  will  emphasize  in  this  paper. 

First,  the  theory  enables  us  to  define  a  real,  symmetric  matrix  Y(E)  that  can 
be  analytically  continued  across  thresholds  and  thereby  relate  the  properties 
of  open  channel  scattering  states  above  threshold  and  predissociating  closed 
channel  resonances  below  threshold  (10,11).  Second,  we  will  explicitly  con¬ 
sider  the  case  of  weak  coupling  between  an  initial  state  and  a  manifold  of 
final  states,  and  show  how  a  half-collision  factorization  of  transition  ampli- 
tudes  for  predissociation  and  photodissociation  can  be  developed  (12).  This 
analysis  will  be  illustrated  by  interpreting  photodissociation  and  radiative 
redistribution  experiments  in  which  an  aligned  atom  is  produced  following 
molecular  absorption  of  a  polarized  photon.  j-j 


2.  MCQDA  FORMULATION 

In  general  the  wavef unction  which  describes  the  diatomic  species  AB  at  total 
energy  E  can  be  expanded  as 

ip^,(E)  -  E  la>  F^^,(E,R)/R  (1) 

a  a  aa 


where  the  basis  functions  la>  span  all  coordinates  except  the  AB  interfragment 
separation  R.  This  may  also  be  written  in  the  following  matrix  notation,  where 
I a>  is  a  row  vector; 

♦  -  la>  F(E,R)/R  (2) 


Introduction  of  eq.  1  into  the  Schrodlnger  equation.  Hip  -  Eip,  generates  the 
standard  matrix  Schrodlnger  equation,  or  equivalently,  the  set  of  close  coupled 
equations , 


2u  . 

^  F  >  CE-I’ 
dR*^ 


W(R)]  F  -  0 


(3) 


where  1®  is  the  unit  matrix  [In  this  paper  diagonal  matrices  are  always  indi¬ 
cated  by  a  superscript  zero].  The  symmetric  normally  nondiagonal  matrix  W, 
which  describes  the  potentials  and  couplings  among  the  various  states  a,  is 
found  by  taking  matrix  elements  of  the  total  Hamiltonian  over  the  basis  func¬ 
tion  la>. 

Although  it  is  possible  to  set  up  and  solve  such  equations  using  a  variety 
of  choices  of  basis  |o>,  the  scattering  boundary  conditions  in  the  limit  of 
large  R  must  be  applied  using  a  particular  basis  set,  namely,  the  channel  state 
basis  |'r>.  These  basis  states  are  eigenstates  of  total  angular  momentum,  ^  « 

J  *  J,  where  J  -  +  Jg  is  total  fragment  angular  momentum  and  t  is  the 

relative  angular  momentum  of  A  and  B.  For  the  case  of  atom-atom  scattering 
this  basis  corresponds  to  Hund's  case  (e)  of  diatomic  angular  momentum 
coupling.  The  potential  matrix  W*  in  the  channel  state  representation  is 
asymptotically  diagonal,  having  elements 


W  (R--) 
YY 


2^ 


(H) 
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The  channel  energy  is  the  total  asymptotic  fragment  energy,  not  Including 
kinetic  energy. 

The  expansion  basis  must  always  be  truncated  with  only  a  finite  number  of 
states,  hopefully  chosen  to  represent  the  physics  adequately.  The  total  number 
of  channels  is  divided  into  Nq  open  and  Nq  closed  channels,  depending  on 
total  energy.  The  channel  wave  vector 

k^  -  C2u(E-E*)/'h2]^/2 

is  real  for  open  channels,  for  which  the  fragments  can  separate  asymptotically 
with  kinetic  energy  E-E^  .  But  k^  is  pure  imaginary  for  closed  channels,  l.e., 

E  <  E^  ,30  that  the  fragments  can  not  separate  in  channel  Y.  If  E  <  E*  for 

di 

all  Y,  -  Nq  and  the  only  possible  states  of  AB  are  bound  states.  If  E  > 
for  all  Y,  N.J  -  Ng  and  inelastic  collisions  can  couple  any  two  channels.  If  E 
is  such  that  some  channels  are  open,  some  closed,  then  it  is  possible  to  have 
mecastable  bound  states,  or  resonances,  which  can  decay  by  predissociation. 

The  MCQDA  gives  a  unified  way  to  describe  open  and  closed  channel  properties 
and  gives  new  insights  into  the  decay  of  predissociating  resonances.  The  basic 
idea  of  the  analysis  is  to  introduce  a  set  of  diagonal  reference  potentials  W® 
such  that 

W(R)  -  ¥®(R)  ♦  V(R)  ,  (6) 

and  7(R)  vanishes  asymptotically.  Each  reference  potential  gives  rise  to 

two  linearly  independent  solutions  of  the  Schrodlnger  equation,  f^  and  g^.  In 
matrix  form  the  diagonal  matrices  f®  and  g®  are  each  solutions  of  the  diagonal 
(or  uncoupled)  set  of  equations 

d^  2u 

— y  F®  *  -y  (E‘1®-¥®)  F®  -  0  .  (7) 

dR^ 

The  general  solution  to  the  coupled  equations  3  can  be  written  as  a  linear 
combination  of  these  reference  solutions,  with  R-dependent  coefficients  A  and 

B, 

F(E,R)  -  ?®(E,R)  A(E,R)  *  g®(E,R)  B(E,R)  .  (3) 

If  ?®  and  g®  are  carefully  chosen  to  be  analytic  in  S,  then  the  N^xN^ 
dimensional  matrix  defined  by  the  asymptotic  quantity 

Y(E)  -  B(E,-)  ACE,-)"’  (9) 

can  be  analytically  continued  across  dissociation  thresholds  (i.e.,  like  a 
"quantum  defect")  and  quantitatively  relates  the  scattering  and  bound  state 
properties  of  t  (S,R).  In  addition,  the  complex  matrix 

t-(E)  -  [A(E,«)±  iB(E, -)]'''  (TO) 

may  be  used  to  define  "hzdf  collision"  amplitudes  which  describe  rinal  state 
distributions  in  photodissociation  and  weak  predissociation  processes.  Note 
that  no  approximations  are  involved  in  writing  eq.  3.  In  fact,  the  original 
second  order  differencial  equations  3  can  be  rewritten  as  two  coupled  ♦’irst 
order  differential  equations  for  A  and  B,  or  equivalently,  for  R-dependent 
analogs  of  Y  and  t“  (9-12).  Thus  our  analysis  is  related  to  amplitude  density 
and  phase  amplitude  methods  in  scattering  t.heory  (20,21).  Our  contribution  is 
to  emphasize  the  intrinsic  analytic  properties  of  the  wavefunction  that  can  be 
projected  from  such  an  analysis. 


i»> 


Although  the  MCQDA  la  a  nonperturbatlve  analysis,  its  usefulness  will  often 
depend  on  making  a  Judicious  choice  of  reference  potentials  W**.  The  reference 
basis  |a>  is  derived  from  the  channel  state  basis  |Y>  by  an  drthorgonal 
transformation 


|o>  -  1T>  i(R) 

H-M*  ♦V-ilT 


(11) 

(12) 


Normally  the  transformation  m  will  be  chosen  for  physical  reasons  to  suggest 
reference  potentials  W**  which  isolate  particular  effects  by  making  V  small  or 
vanishing  at  some  value  of  R.  For  example,  if  we  want  to  be  able  to  generate 
an  energy  Insensitive  Y(E)  for  extrapolating  across  thresholds,  we  would  seek  a 
transformation  that  gives  a  coupling  V(R)  that  operates  only  in  a  short  range 
zone  of  R,  typically  in  the  R^  region  of  Born-Oppenheimer  potentials.  On  the 
other  hand,  for  describing  a  photodissociation  process,  we  would  want  to  choose 
a  V  which  is  negligible  in  the  zone  of  R  near  Rg  which  determines  the 
Franck-Condon  absorption,  whereas  there  may  be  strong  inelastic  couplings  due 
to  V  in  an  intermediate  range  of  R  that  determine  the  distribution  of  atomic 
fragment  states.  Depending  on  the  specific  application,  some  choices  of 
suitable  reference  functions  for  diatomic  processes  could  be:  a  pure  Hund's 
coupling  case,  such  as  the  standard  case  (a)  states;  a  set  of  intermediate 
coupled  states;  an  adiabatic-electronic-rotational  (AER)  basis,  which 
continuously  diagonalizes  the  electronic-rotational  Hamiltonain  as  a  function 
of  R  and  smoothly  transforming  from  one  Hiuid’s  case  to  another  (8b, 9). 

Several  types  of  reference  solutions  may  be  defined,  depending  on  boundary 
conditions  chosen  (10,18).  The  solutions  useful  for  extrapolating  across 
thresholds  are  the  functions  ?  and  g,  respectively  regular  and  irregular  as  R  + 
0.  They  are  defined  quantum  mechanically  but  have  the  following  semiclassical 
interpretation  in  a  classical  region  of  phase  space: 


sin  Sg(R)/ky^(R)  and 


cos  85,(R)/ky^(R) 


(13) 


where  S~(R)  is  the  WKB  phase  angle  (a  is  the  classical  turning  point) 


(14) 


Each  reference  potential  ^^^^(R)  uniquely  defines  two  energy- dependent  phase 


angles  ^^^(E) 


and  Vg(E) 


Above  a  channel  threshold  the  open  channel  solutions 
define  the  reference  phase  shifts  (we  ignore  subtle  threshold  effects  discussed 
in  detail  in  references  10  and  11) 


(R— )  -  sin(k^R-  f  +  C,,) 


(15) 


Below  a  channel  threshold  the  regular  solutions  fjj(E.R)  are  asymptotically 
divergent  except  at  closed  channel  eigenvalues  E  »  E„  defined  by  the  modular  it 


values  of  Vjj(E) 


tan  V  ( E)  "0  or  v_(E*)  ■  nir 
a  on 


(16) 
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The  semiclassical  approximation  to  v 


Va  ♦  Y  -  Z**  kg(R)  dR 


(17) 
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gives  the  standard  Bohr-Soomerfeld  quantization  rules. 

The  MCDWA  amalysis  of  Mies  (10)  shows  how  to  define  the  symmetric  real 
matrix  T(E)  in  eq.  9.  If  all  channels  are  open,  the  N.pXN>p  scattering  matrix  S 
is  simply  related  to  Y  as  follows 

S  -  (1»  ♦  iY)  (1»  -  il)"''  e^5“  .  (i=S) 

If  we  calculate  S  and  the  reference  phase  shift  diagonal  matrix  5®  above 
threshold,  then  we  can  calculate  the  N.[.xN-r  Y  matrix.  If  Y  is  sufficiently 
energy  insensitive,  then  it  may  be  extrapolated  below  thresholds  to  predict  the 
resonance  scattering  and  closed  channel  properties.  For  example,  if  we  take 
the  case  N<p-2  with  channel  1  closed  and  channel  0  open,  then  the  elastic 
scattering  phase  shift  of  the  open  channel  is  modified  due  to  predissociation 
originating  from  the  bound  states  located  at  v^(E)  «  nir  (11), 

Y  Y 

no(E)  -  £o(E)  +  tan’’  CYqo  -  - ]  (19) 

CHfi  (E)“Y-|  "I 

If  both  channels  are  closed,  the  bound  state  eigenvalues  are  given  by  the 
expression 

(tan  YQ(En)  +  Yqq)  (tan  ^^(En)  +  Y^,)  -  Y^^  (20) 


Eqs.  19  and  20  are  valid  for  arbitrarily  strong  coupling  and  are  only  limited 
by  our  ability  to  obtain  the  reference  phases  and  the  matrix  Y(E)  -  constant. 

In  the  limit  of  weak  coupling  when  Y  ♦  0,  eq.  (19)  transforms  into  the 
usual  Lorentzian  expression  for  an  Isolated  resonance.  The  weak  coupling  width 
Fjj  and  shift  AEj^  of  the  predlssociatlng  level  at  E  -  Ejj  are 

1  2AGn  2  1  AG_ 

I-,  - - a  Y  and  AE  - - Q-Y,,  ,  (21) 

n  ir  01  n  ir 

where  we  have  Introduced  the  approximation 


3E  E-E^ 


(22) 


where  AGjj  is  the  vibrational  spacing. 

We  have  carried  out  a  number  of  test  calculations  on  2  state  models,  i.s.,  2 
open  or  1  open  and  1  closed  channels,  and  have  demonstrated  the  success  of  the 
analysis  over  a  wide  range  of  conditions  (11).  These  include  an  analysis  of 
Cd2  state  predissociation  to  Cd(3pQ)  Cd(’SQ)  using  AER  reference  states, 
and  an  analysis  of  the  O2  '^1^  predissociation  using  diabatic  reference  states. 
Model  strong  coupling  cases  based  on  the  latter  were  used  to  demonstrate  the 
nonperturbative  nature  of  the  analysis.  The  Y(E)  matrix  was  found  to  be  most 
energy  insensitive  for  short  range  inner  potential  wall  curve  crossings  near 
crossings  where  two  reference  potentials  remain  nearly  parallel  over  a  r..-ige  of 
R.  Its  greatest  energy  sensitivity  will  occur  for  sharp  outer  crossings  due  to 
the  more  rapid  energy  variation  of  Franck-Condon  overlaps. 


3.  WEAK  COUPLING  HALF  COLLISION  ANALYSIS 


Let  us  now  consider  a  special  case  which  covers  a  wide  variety  of  actual 
phenomena  and  for  which  the  MCQDA  allows  us  to  develop  a  number  of  useful 
insights  or  approximations.  Assume  that  we  have  a  manifold  of  Nj  channels 


which  may  be  strongly  coupled  among  themselves  and  a  similar  manifold  "of  Np 
final  channels,  but  that  the  coupling  between  the  Initial  and  final  manifolds, 
described  by  the  NpXNj  Vpi  submatrix.  Is  very  weak.  For  simplicity  we  will 
consider  the  case  here  where  Nj  ■  1 ,  l.e.,  there  Is  only  a  single  nondegenerate 
Initial  channel.  “If  the  Initial  channel  Is  open,  as  well  as  the  final  chan--  ' 
nels,  the  free-free  transition  amplitudes  are  given  by  the  NpXNj  subblock  of 
the  S  matrix,  described  to  first  order  In  Vpj  by  the  Born  distorted  wave 
approximation, 

Spj(E)  -  -2irl  <*f(E)|Vl*J(E)>  .  (23) 

when  tp*  are  energy  normalized  wavefunctlons  calculated  In  the  absence  of  the 
weak  coupling  and  the  ±  refer  to  incoming  and  outgoing  scattering  boundary 
conditions.  Similarly,  If  the  Initial  channel  la  closed,  the  bound-free 
predlsaoclatlon  transition  amplitudes  for  decay  of  the  closed  channel  resonance 
at  E-Ejj  are 


TpiCEn) 


-2irl  <*^(En)|V|En> 


where  lEjj>  Is  the  unit  normalized  bound  state  wavefunction. 

There  are  numerous  examples  of  weak  predissociations  In  molecular  spectros¬ 
copy.  One  special  case  Is  predissociation  induced  by  a  radiation  field, 
namely,  photodissociation.  If  the  initial  state  basis  functions  are  molecule- 
radiatlon-field  states  | aj>  |nqu>,  with  n  photon  of  frequency  ui  and  polariza¬ 
tion  q,  and  the  final  states  have  n-1  photons  |op>ln-1 ,q,to>,  then  the  free-free 
process,  eq.  23,  describes  line  broadening  and  radiative  redistribution 
phenomena  and  the  bound-free  process,  eq.  24,  describes  photodissociation  (22). 
For  the  case  of  radiative  coupling  the  weak  coupling  matrix  is  Just  (22,23). 


■(rad) 


-]  <aF|eq*ulai;> 


where  ^  is  the  field  Intensity  (photons  cm'^sec"'),  e^  the  polarization 
vector  of  the  light,  and  u  the  molecular  dipole  operator. 

When  we  impose  scattering  boundary  conditions  on  the  asymptotic 
wavefunction  in  the  channel  state  basis  |T>,  we  find  (using  eq.  11)  that  the 
energy  normalized  wavefunctlons  for  the  respective  initial  and  final  manifolds 


wnere  k  -  (2u/fi^ir)^‘^^  and  the  normalization  matrix  is 

N±  -  Ct*(-)]"’  ,  (27) 

with  t  given  by  eq.  10.  (12)  We  now  wish  to  make  one  more  assumption.  Assume 
that  the  weak  coupling  Vpj  which  couples  the  two  manifolds  is  only  significant 
in  some  short  range  zone  of  R  near  R^  and  that  the  final  state  interactions  Vpp 
which  determine  the  final  state  distribution  are  important  in  a  longer  range 
zone  of  R.  Since  B  must  vanish  at  the  origin  and  does  not  become  large  until 
R  >  Rg,  the  weak  transition  occurs  in  a  zone  of  R  where  to  a  good  approximation 
we  may  take  A-1,  B-0.  We  immediately  see  from  eqs.  23,  24,  26  that  the  weak 
coupling  amplitudes  are 

Spj  -  -2Tri  M^(E)  Spi(E)  N*(E)  (28a) 


(28b) 


where 
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and  we  have  used  the  result  -  N“*. 

The  results,  eq.  28,  show  how  a  very  useful  factorization  of  the  weak  cou¬ 
pling  transition  amplitudes  can  be  achieved.  The  3  and  t  amplitudes  are  real 
Franck-Condon  amplitudes  which  describe  the  excitation  process.  The  Mp  matrix 
is  a  "half  collision"  dynamical  matrix  which  incorporates  the  effect  of  final 
state  interactions  and  provides  the  correlation  between  the  short  range 
reference  states  excited  by  Vpj  and  the  asymptotic  distribution  of  fragments 
which  could  be  measured  experimentally.  The  initial  state  matrix  Mj  is  a  sim¬ 
ple  phase  factor  for  a  single  initial  channel,  whereas  for  a  manifold  of  ini¬ 
tial  channels,  it  would  describe  the  formation  of  the  reference  states  from  the 
asymptotic  states.  We  see  that  the  physical  picture  implicit  in  the  factor¬ 
ization,  eq.  28b,  is  the  time-independent  analog  of  the  time-dependent  wave- 
packet  viewpoint  of  photodissociation  of  Lee  and  Heller  (24).  According  to 
this  viewpoint,  photodissociation  is  described  as  the  initial  creation  of  an 
excited  state  wavepacket  through  a  "Franck-Condon"  excitation  from  the  initial 
bound  state  due  to  weak  radiative  coupling.  The  time  evolution  of  this  wave- 
packet  leads  to  the  distribution  of  final  states. 

We  have  shown  two  Important  properties  of  the  matrix  (12).  First,  as 
long  as  no  nonclassical  regions  are  encountered  by  the  fragments  between  the 
excition  region,  R  »  R^,  and  the  detection  region,  R  -  then  H'*'  is  unitary  to 
a  very  good  approximation.  Second,  the  matrix  can  be  further  factored  as 

-  0^  O'"  (30) 

The  factor  U*,  the  so-called  frame  transformation  matrix,  depends  only  on  the 
reference  solutions  and  the  orthogonal  basis  transformation  matrix  a.  If  an 
AER  reference  basis  is  used,  it  takes  on  the  form 

-  L(Rg)  e^^’  (31) 

with  L(R)  -  e^®*"  m(R)  e"^®"  (32) 

where  8®  is  the  matrix  of  WK3  reference  phases,  eq.  14.  The  distance  Rg  is  a 
characteristic  switching  distance  at  which  the  final  state  interactions  begin 
to  become  important  as  the  fragments  separate.  The  matrix  0"^  in  eq.  30  incor¬ 
porates  the  effect  of  final  state  interactions  in  causing  deviations  from  the 
pure  frame  transformation,  eq.  31 •  It  is  possible  to  obtain  a  first  order 
differential  equation  for  0*  by  which  these  deviations  can  be  calculated  (12). 

A  special  limiting  case  which  has  been  used  to  predict  final  state 
distributions  in  diatomic  photodissociation  is  the  recoil  approximation,  or 
more  specifically,  the  J-independent  recoil  approximation  (25).  This  amounts 
to  ignoring  any  J  dependence  of  final  state  Franck-Condon  factors  and  phase 
development.  This  limit  will  always  apply  at  sufficiently  high  fragment 
kinetic  energy.  In  our  formulation,  this  limit  corresponds  to  setting  * 
m(RQ).  If,  however,  we  set  -  U'^,  thereby  ignoring  the  nonadiabatic 
corrections  in  0'^  but  incorporating  the  effect  of  adiabatic  phase  development, 
and  continue  to  ignore  any  J-dependence  in  the  Franck-Condon  amplitudes,  §  or 
t,  we  obtain  an  improved  approximation  which  allows  us  to  incorporate  the 
effect  of  molecular  axis  rotation  into  orientation/alignment  experiments.  We 
thus  see  that  the  factorizations,  eqs.  28  and  30,  are  very  powerful  tools  which 
enable  us  to  introduce  a  heirarchy  of  approximations  for  successively  introduc¬ 
ing  various  dynamical  effects  into  the  interpretation  of  photodissociation  or 
predissociation  product  distributions. 
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■  In  our  earlier  discussion  we  emphasized  the  Y  matrlx"and  Its  continuity - 

across  a  dissociation  threshold.  For  our  present^example  of  a  single  initial 
channel  weakly  coupled  to  a  manifold  of  final  channels,  the  Ypj  submatrix  to 
first  order  in  Vpj  is 


Ypj(E)  -  -ir  A^\e,»)  SpiCE)  . 


(33) 


The  continuity  of  Y  across  the  initial  state  threshold  Is  thus  reflected  by  the 
well  known  continuity  of  Franck-Condon  amplitudes  The  connection  between 
and  tpj  In  eq.  29  Is  readily  seen  when  we  write  the  unit  normalized  Initial 
bound  state  function  as 


lEn>  -  ^  (34) 

where  we  have  used  eq.  22  and  the  results  of  references  3  and  10.  Thus, 

tpi(E„)  -  (j^)  (35) 

shows  the  relation  between  the  nonanalytic  t  and  the  analytic  §. 

One  important  consequence  of  the  unitarity  of  M*  is  that  the  total  predis¬ 
sociation  decay  rate,  given  by  the  following  scalar  dot  product, 

Aj^  ”  Tp  Tp/h  •  (4ir^  ^FX  ^FI^  (AGjj/h),  (36) 

depends  only  on  the  excitation  Franck-Condon  factors  §,  and  not  on  the  half¬ 
collision  dynamical  matrix.  This  justifies  the  usual  adiabatic  approximation 
for  continuum  spectroscopy  and  weak  predissociation,  i.e.,  given  our  assump¬ 
tions,  it  is  only  necessary  to  take  into  account  the  local  interactions  in  the 
Franck-Condon  region  to  calculate  the  total  rate  of  decay.  Since  when  the 
initial  channel  is  open  the  first  factor  in  parenthesis  in  eq.  36  is  inter¬ 
preted  as  the  total  inelastic  transition  probability  from  I  to  the  F  manifold, 
we  see  that  the  closed  channel  decay  rate  is  simply  interpreted  as  a  single 
cycle  transition  probability  times  the  oscillation  frequency  AG^j/h.  The 
continuity  of  §  across  threshold  and  this  interpretation  of  eq.  36  has  recently 
been  discussed  by  Tellinghuisen  and  Julienne  (26)  in  relation  to  the  inverse  of 
the  phenomena  considered  here,  namely,  for  radiative  association  and  bound 
state  radiative  decay. 

We  have  carried  out  specific  calculations  for  the  case  where  polarized  ab¬ 
sorption  from  a  ^Z'*'  state  through  either  a  ^  Z  or  ^  II  molecular  Born-Oppenheimer 
reference  state  leads  to  fragment  atoms  (23,27).  The  ^P  atom  may  be 

aligned  (q-0)  or  oriented  (q»±1).  This  describes  both  Sr  or  Ba  +  rare  gas 
radiative  redistribution  experiments  (28)  and  the  Ca2  photodissociation  experi¬ 
ment  of  Vigu4,  et.  al.  (29).  Our  full  close  coupled  results  for  Sr  +  Ar  are 
in  good  agreement  with  the  radiative  redistribution  experiments.  The  half  col¬ 
lision  analysis  leads  to  an  especially  simple  and  pleasing  physical  picture  of 
the  product  alignment. 

The  same  angular  momentum  transfer  formalism  applies  to  both  radiative 
redistribution  (23)  and  photodissociation  (30).  The  transfer  cross  sections 
require  a  coherent  sum  of  transition  amplitudes  resulting  from  P,  Q,  and  R 
type  radiative  transitions,  i.e,  those  having  AJ  »  1,  0,  and  +1  respectively. 
The  transfer  cross  sections  that  result  from  applying  the  factored  form  of  the 

matrix  can  be  expressed  in  terms  of  an  angle 

-  Ag/2  . 
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The  angle  aj^/2  originates  from  the  frame  transformation  matrix  U'*'  in  eq.  30  and 
represents  the  effect  of  adiabatic  motion  on  the  Born-Oppenheimer  potential  A. 
Semiclassically,  it  is  the  angle  through  which  the  molecular  axis  rotates  in  a 
complete  half  collision  in  the  final  state,  that  is,  (x-0q)/2.  where  is  the 
final  state  classical  deflection  function.  The  classical  deflection  function  - 
enters  since  it  is  23^y^/3J  and  the  adiabatic  phases  vary  slightly  with  Jp 
for  adiabatic  evolution  following  P,  Q,  and  R  excitation.  The  angle  A3/2 
originates  from  the  effect  of  nonadiabatlc  Coriolis  coupling  between  I  and  II 
states  in  the  nonadiabatlc  correction  matrtix  0*^.  It  represents  the  decoupling 
of  electronic  amgular  momentum  from  the  rotating  axis  beyond  a  switching  dis¬ 
tance  Rg. 

This  model  leads  to  the  same  polarization  ratio  as  found  by  Lewis,  et  al. 
(31)  on  the  basis  of  an  ad  hoc  geometric  model.  This  model  corresponds  to  the 
picture  that  angular  momentum  projection  A  rotates  with  the  Internuclear  axis 
inside  of  some  characteristic  Internuclear  separation  Rg  but  is  decoupled  and 
fixed  in  space  for  R  >  Rg.  If  molecular  axis  rotation  is  neglected,  our  half 
collision  analysis  predicts  the  j-independent  recoil  limit  results  (25,29), 
namely,  respective  q  -  0  polarization  ratios  of  1/2  and  ?/9  for  Z  and  H  exci¬ 
tation.  The  usefulness  of  the  complete  half  collision  analysis  is  that  correc¬ 
tions  to  this  simple  J-lndependent  recoil  approximation  are  possible  by  intro¬ 
ducing  additional  effects  into  Np.  In  the  present  case  it  happens  to  be  simple 
to  isolate  the  separate  effects  of  adiabatic  and  nonadiabatlc  dynamics. 

We  have  also  used  the  half  collision  emalysis  to  interpret  Na  *■  rare  gas 
radiative  redistribution  experiments  (31)  in  which  the  Na  ^P(j-1/2)/^P(j-3/2) 
ratio  was  measured  following  wing  absorption  to  the  and  states  (32).  The 
recoil  limit  result  is  approached  for  fragment  kinetic  energy  of  about  2000 
cm“^  for  an  argon  collision  partner,  but  at  a  much  lower  kinetic  energy  for  the 
lighter  He  collision  partner.  The  large  departures  of  the  fine  structure 
branching  ratio  from  the  recoil  limit  for  low  (thermal)  separation  energies 
reflects  the  influence  of  the  adiabatic  potentials  and  their  curve  crossings. 
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ABSTRACT 


For  dipole  allowed  atomic  transitions  the  radiative  matrix  element 
which  defines  the  pressure-brcadened  atomic  llneshape  is  only 
conditionally  convergent.  Using  a  commutator  technique  to  redefine  the 
Integral,  we  isolate,  and  ultimately  reject,  the  contribution  of  an 
indeterminant  asymptotic  surface  integral  that  is  associated  with  the 
energy  normalization  of  the  continuum  wavefunctlons  which  describe  the 
binary  collision  of  the  atom  and  its  perturber.  The  remaining 
contributions  which  are  absolutely  convergent,  give  the  multichannel 
atomic  llneshape  which  includes  effects  of  non-adlabatlc  and  inelastic 
scattering.  Further  we  show  the  relationship  of  the  commutator  integral 
to  the  exact  requirements  of  close-coupled  scattering  theory  for 
radlatively-induced  collisions.  This  scattering  analysis  suggests  the 
interpretation  of  the  convergent  llneshape  as  an  expression  of 
multichannel  Inelastic  collisions  between  field-dressed  atomic  states. 
This  same  interpretation  applies  both  in  the  impact  and  the  static  limit. 
In  this  paper  we  make  explicit  comparisons  which  demonstrate  the 
equivalence  between  the  commutator  integral  and  the  numerical 
close-coupled  results  in  the  weak-^ield  limit.  Here  we  emphasize  the 
static  limit,  well  in  the  wings  of  the  atomic  line,  where  the  Jablonski 
stationary  phase  WK3  analysis  is  often  applied  with  good  effect. 


1 .  Introduction 

The  calculation  of  the  pressure-broadened  lineshape  for  the  dipole 
allowed  atomic  transition  hu  Af*  requires  the  evaluation  of  the 
following  continuum-continuum  matrix  element  (Mies  1981;  Julienne  and  Mies 
1984,  Van  Regemorter  and  Feautrier  1985) 

-  -2x1  <f“(E")ll(a.)*il'l'*(E*)>  .  (1) 

where  the  field-free  multichannel  wavefunctions  f-  will  be  defined  in 
section  2  below.  This  expression  applies  in  the  weak-fleld  limit,  when  the 
electric  field  strength  E((o)  associated  with  the  photon  hu  is  very  small 
and  the  radiative  scattering  matrix  element  can  be  evaluated  in 

perturbation  theory.  The  precise  details  of  the  lineshape  analysis  depends 
on  the  extent  to  which  the  incident  photon  is  detuned  from  the  exact  atomic 
resonance  energy  hu^^  -(E^  “El),  where  the  .detuning  parameter  Afi  is, 

o  * 

Vi  •  (hai-hu»^j^)  ■  (hu-Ep  El)  .  (2) 

In  the  wings  of  the  line,  where  is  very  large  compared  to  the 
homogeneous  width  of  the  atomic  transition,  the  absorption  coefficient  is 
directly  proportional  to  the  radiative  scattering  cross  section  determined 
by  the  transition  matrix  elements  in  Eq.  (1)  (Julienne  and  Mies  1984). 

The  Lorentzian  linewidth  given  by  the  impact  approximation  in  the  vicinity 

part 

of  the  line  center  can  also  be  expressed  in  terms  of  3^.,  ^Mies  and 
Julienne  1985,  Van  Regemortor  and  Feantrier  1985). 


The  expression  of  the  llneshape  In  terms  of  Eq.  (1)  is  a 
generalization  of  the  standard  time- independent  molecular  Franck-Condon 
factor  approach  to  the  line  profile.  A  good  recent  review  of  line 
broadening  theory  which  discusses  the  interrelationships  and  historical 
development  of  the  various  techniques  has  been  given  by  Allard  and 
Kielfopf  (1982).  The  numerical  evaluation  of  the  Integral  (1)  for  dipole 
allowed  transitions  requires  some  care  since  the  transition  dipole 
operator  fi  is  non-vanishing  even  when  the  absorbing  atom  A  and  the 
perturber  B  are  infinitely  separated.  Thus  the  resultant  oscillatory 
integrand  in  Eq.  (1)  persists  forever  and  does  not  give  an  absolutely 
convergent  Integral,  and  suitable  cechnlques  must  be  used  to  extract  the 
physically  meaningful,  absolutely  convergent  components  in  this  matrix 
element. 

The  standard  way  aro'ind  this  problem  is  to  use  Jablonskl’s  method 
(1948).  By  invoking  the  stationary  phase  approximation  with  JWKB  wave- 
f'lnctlons,  the  offending  terms  are  eliminated  by  a  random  phase  argument, 
and  simple  semlclasslcal  approximations  to  Eq.  (1)  result.  These  provide 
the  basis  ^or  the  unified  Franck-Condon  theory  of  the  profile  (Szudy  and 
Bayllss  1975)  and  also  provide  Justification  for  the  standard  quasistatic 
expressions  for  the  wing  profile  (Hedges  et  al  1972;  see  also  Herman  and 
Sando  1978). 

Convergent  niimerlcal  techniques  have  been  devised  for  the  direct 
numerical  evaluation  of  the  full  quantum  integral  (1).  This  has  been  done 
in  two  Independent  ways  for  neutral  broadening  by  Herman  and  Sando  (1978) 
and  by  Bleniek  (1981).  However,  their  fornwlae  have  been  restricted  to  the 


adiabatic  approximation,  that  Is,  to  using  single  channel  wavef<jnctlons  In 
(1).  Actually  convergent  ways  to  evaluate  (1)  for  multichannel 
wavefunctlons  have  already  been  derived  for  application  to  electron-atom 
(Feautrler  et  al  1976)  and  electron-ion  (Seaton  1984)  broadening  problems. 

In  this  paper  we  show  how  the  commutator  method  of  Herman  and  Sando 
(1978)  can  be  generalized  to  yield  an  absolutely  convergent  Integral  using 
multichannel  wave f'inct ions.  The  discussion  is  given  in  terms  of  the 
simple  case  of  broadening  of  an  atomic  transition  by  collllslon  with  an 
atomic  perturber.  We  emphasize  that  our  procedure  is  applicable  to  more 
complicated  systems.  For  example,  our  formalism  can  be  generalized  by  a 
reinterpretation  of  channel  labels  to  treat  broadening  of  molecular 
adsorption  lines  by  perturbers  with  Internal  stnicture.  Our  general 
expression  reduces  to  those  of  Herman  and  Sando  (1978)  and  Bieniek  (1981) 
as  special  cases  for  single  channel  wavefiinctlons.  We  also  show  how  a 
distorted  wave  analysis  of  the  radiative  scattering  matrix  for  dressed 
atoms  ( Cohen-Tannoudj 1  and  Reynaud  1977,  Courtens  and  Szbke  1977,  Mies  and 
Ben  Aryeh  1981,  Mies  1984)  lead  to  the  same  formula.  Thus,  either 
close-coupled  scattering  calculations  in  which  the  radiation  field  couples 
asymptotically  well-defined  dressed  states,  or  numerical  evaluation  of  the 
perturbative  expression  (1)  can  be  used  to  calculate  the  nonadlabatic 
version  of  the  molecular  Franck-Condon  Integrals. 

The  next  section  introduces  the  notation  and,  together  with  Appendix 
A,  defines  the  multichannel  wavefunctlons.  The  third  and  fo'irth  sections 
respectively  describe  the  commutator  and  dressed  atom  scattering 
derivations  of  convergent  Integral  fonvilae  which  replace  (1).  The  final 


section  gives  some  Illustrative  numerical  calculations  and  a  brief 
discussion  of  llneshape  theory  in  the  wings  of  pressiire  broadened  lines 


The  core  profile,  with  special  emphasis  on  the  impact  limit,  together  with 
some  consideration  of  strong-field  effects  will  be  addressed  in  an 
accompanying  paper  (Mies  and  Julienne  1985). 

2.  Single  perturber  Theory 

We  consider  broadening  of  the  atomic  transition  line  by  a 
perturber  B  in  the  binary  collision  limit.  In  the  so  called  optical 
collision  process  (Mies  1981)  of  interest  here,  the  tr2msition  is 
asymptotically  allowed.  The  two  collision  partners  are  treated  as  one 
absorbing  diatomic  species. 


+  *  Il(i»  CA*+*»B))f  ♦  Ef] 

(3) 

>  floj  E" 

(4) 

The  labels  "i"  and  "f"  denote  the  internal  state  of  the  quasidlatomlc  pair 
(with  asymptotic  internal  energy  and  Ef)  before  and  after  the  binary 
absorption  event,  respectively.  The  superscript  refers  to  the  system 
of  incoming  collision  partners  with  total  energy  E"*”,  while  indicates 
outgoing  particles  after  absorption  of  a  photon  during  a  collision.  The 
multichannel  wavefunctlons  1'^(E'*')  and  ^^(E")  that  appear  in  the  matrix 
element  (1)  describe  field-free  continuum  motion  of  A  and  A*, 
respectively,  with  respect  to  3.  The  wave  f'inctlon  'F^(E'*')  contains 
incoming  particle  flux  in  state  i  and  outgoing  particle  flux  in  all 
channels.  1'T(E~)  has  incoming  particle  flux  in  all  channels,  but  outgoing 


flux  In  state  f  only.  [The  properties  of  multichannel  continuum  wave 
functions  have  been  extensively  discussed  by  Mies  (1980,  1981,  1984)  and 
Singer  et.  al.  (1982.  1984.  1985).] 

For  simplicity  of  presentation  we  consider  the  collision  partner  B  to 
be  structureless.  In  this  case,  Ej(Ef)  is  the  internal  energy  of  the 
A(A*)  fragment.  Generalization  to  more  complicated  cases  in  which  the 
internal  state  of  either  atom  may  change  during  an  optical  collision  is 
accomplished  by  making  the  labels  i  and  f  composite  labels  that  indicate 
the  asymptotic  energy  levels  of  both  fragments.  The  same  is  true  for 
optical  collisions  involving  molecular  species,  in  vrtiich  the  channel 
labels  differentiate  rovlbronic  levels.  In  order  not  to  obscure  the 
presentation  of  oiir  method,  we  treat  a  very  simple  physical  system  in  our 
derivations. 

As  Che  interatomic  radial  coordinate  R  becomes  infinite,  the  fragments 
are  further  characterized  by  an  Incoming  e*  and  outgoing  asymptotic 
relative  kinetic  energy  which  defines  the  total  energies  associated  with 
the  Incoming 

e*  (5a) 

and  outgoing 

E  ■  ••  E^  ♦  Au  ,  (5b). 


molecular  wavef'onctions  respectively. 

The  conservation  of  energy  dictated  by  Eqs.  (4)  and  (5b)  indicates  chat  an 
appropriate  change  in  kinetic  energy  compensates  for  any  detuning  of  the 
absorbed  photon. 
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The  Internal  energy  levies  of  the  A***B  or  quatsimolecular  pair  can 

be  described  by  a  variety  of  equally  acceptable  R-dependent 
electronic-rotational  basis  functions.  However,  to  describe  proper 
scattering  boundary  conditions  It  Is  Important  to  choose  basis  functions 
which  asymptotically  correlate  with  well-defined  atomic  levels  of  the 
non-lnteractlng  fragments.  Such  basis  functions  are  referred  to  as  channel 
states. 

In  general,  the  A-^*  transition  in  Eq.  (3)  will  Involve  a  set  of 
1»1,N'*‘  Initial  and  a  separate  set  of  f»1  ,N“  final  molecular  channel 
states  which  correlate  with  the  product  of  atomic  states  associated  with 
the  A***B  and  A*«*B  collision  complex  respectively.  We  group  these  sets 
Into  two  row  vectors 

4*"  ■  i>J,****.  (7a) 

,*•••.  <P  ,••••,  <*  _]  (7b) 

1  f  N 

These  channel  states  are  but  finite  sub-sets  of  f'inctions  which  span  the 
entire  space  of  the  molecular  Hamiltonian  defined  by  the  collision 
complex,  with  the  single  exception  of  the  radial  coordinate  R.  They  are 
constructed  with  R-dependent  molecular  electronic-rotational  wavef'inctlons 
which  are  combined  to  insure  that,  asymptotically,  (R-«)  and  (R-«) 
become  exact  eigenfunctions  of  the  atomic  Hamiltonians  (H^  *  Hg) , 
including  any  pertinent  relativistic  operators.  F'orher,  each  initial  and 
final  asymptotic  channel  state  is  an  exact  eigenfunction  of  the 


I 


I 
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relativistic  angular  moinentum  operator  and  If  respectively.  In  the 
absence  of  radiation,  the  total  angular  moaienta  of  the  collision  complex 
are  conserved  and  define  exact  quantum  numbers  M'*’  and  J~,M~  at  all 
Interatomic  distances.  In  the  notation  of  molecular  spectroscopy  (Mies 
1968,  1980,  1984)  such  asymptotic  channel  states  correspond  to  Hund's  case 
(e)  electronic-rotational  states.  They  have  also  been  designated  as  the 
space-fixed  atomic  basis  In  the  half-collision  analysis  of 
photodlssoclatlon  of  Singer,  Freed  and  Band  (1983,  1985). 

Close  coupled  scattering  theory  Is  predicated  on  the  assumption  that 
a  finite  set  of  such  channels  can  be  preselected  which  are  sufficient  to 
describe  any  inelastic  scatterings  involved  In  the  collision.  Two  such 
sets  are  presumed  to  have  been  chosen  In  Eq.  (7),  each  sufficient  to 
describe  the  inelastic  scatterings  of  atom  A  and  A*  respectively.  We  use 
the  first  row  vector  to  expand  the  complete  molecular  wavefunction  at 
total  energy  E* 

f*  (E*)  -  /  f""  (E*,R)/R  (8a) 

and  obtain  a  comparable  row  vector  of  orthogonal,  energy-normalized 
continuum  wavef'inctions  T'*’(E^)  whose  total  degeneracy  N"^  mimics  the  momber 
of  open  (energetically  accessible)  atomic  states  which  are  contained  in 

Quite  independently,  we  obtain  a  second  set  of  continmom  wave- 
f'onctlons  which  we  associate  with  the  final  state  vector 

f“(E')  -  ♦"F"(E“,a)*/R  .  (8b) 

However  the  total  energy  E~  is  prescribed  by  the  incident  molecular  energy 
E^  and  the  incident  photon  in  Eq.  (5b).  Please  observe  that  the 
notation  +  and  -  is  chosen  here  to  Identify  the  separate  seta  of  incoming 
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and  outgoing  continuum  states  Involved  In  the  photon  absorption  process, 
as  well  as  to  label  the  omtually  orthogonal  Initial  and  final  sets  of 
channel  states  which  are  used  to  expand  these  multichannel  wavef’onctlons . 
Actually  both  F'*'  and  F~  are  chosen  to  represent  matrices  of  Incoming 
radial  solution  vectors.  The  final  state  In  Eq.  (1)  Is  required  to 
exhibit  outgoing  boundary  conditions  and  Is  simply  given  by  the  complex 
conjugate  F“(E“,R)*  In  Eq.  (8b).  These  features  are  discussed  In  the 
Appendix.  The  relation  of  to  outgoing  particle  flux  and  final  state 
distributions  In  photodissociation  has  been  discussed  by  Band  e^  ^ 
(1981).  Some  of  the  confusing,  and  conflicting,  aspects  of  the  notations 
used  In  scattering  theory  are  clarified  by  Roman  (1965). 

For  use  In  Eq.  (1),  we  are  required  to  Impose  the  following 
normalization  condition,  independently,  on  each  of  the  square  arrays  of 
radial  solution  vectors  In  Eq.  (8) 

09 

/  dR  f*(E,R)**F±(E=*=,R)  -  6(E“E±)  1*  (9) 

o 


where  1-  denotes  a  N*xM-  dimensional  unit  matrix.  Each  vertical  vector  in 
the  matrix  F*  represents  a  particular  incoming  solution  to  the  appropriate 
set  of  coupled  radial  equations, 

[T^  +  V'‘=(R)  -  (E*,R)  -  0  (10) 
with  asymptotic,  multichannel,  scattering  boundary  conditions  as 


prescribed  in  the  Appendix  A.  Both  the  radial  kinetic  energy  operators 


and  the  matrices  of  asymptotic  relative  kinetic  energies  and 

e“  -  associated  with  the  channel  vectors  and  4"  in  Eq.  (7)  are 

perfectly  diagonal.  Only  the  R-dependent,  real,  symmetric,  interaction 
matrices  Y^(R)  are  non-diagonal,  and,  to  satisfy  our  channel  state 
requirements,  even  these  must  he  constrained  to  vanish  identically  at 
infinite  interatomic  separation. 
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?^(R)  -  9±(R)  -  T±(R)*  -  0  (12) 

R— 

If  radial  couplings  due  to  the  implicit  R-dependence  of  >i(R)  are 
significant,  as  for  instance  when  channels  are  constnicted  with  adiabatic 
electronic-rotational  states  (Mies  1980),  then  antisymmetric  radial 
operator  3/3R  terms  must  be  included  in  ¥^(R) ,  but  the  asymptotic 
constraint  still  persists. 

In  general  there  will  be  large  off-diagonal  elements  in  the 
interaction  matrix  Y*  which  can  be  expected  to  introduce  strong  inelastic 
couplings  among  the  set  of  N*  incoming  channels  Similarly,  the  7~ 

matrix  can  cause  significant  non-adlabatic  mixings  among  the  N~  outgoing 
channel  states  These  inelastic  scattering  effects  are  summarized  by 

the  pair  of  unitary,  symmetric  scattering  matrices  S*(E*)  and  S“(E~)  which 
we  obtain  by  independently  solving  the  two  sets  of  coupled  equations  in 
Eq.  (10)  with  the  bo'indary  conditions  in  Appendix  A  yielding  the  following 
constraints 
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;  dR  #‘^(E,R)  ‘F^CE^.R)  -  5(E-E‘*’)  S'^CE*) 


/  dR  f"(E,R)**F~(E",R)*  -  6(E-E")  S"(E“)' 


For  optical  collisions  the  two  sets  of  channel  states  selected  In  Eq. 
(7)  must  be  of  opposite  total  parity,  and,  thus,  can  never  be  coupled  by 
any  terms  In  the  Isolated  molecular  Hamiltonian  Hj^g(R}, 

-  0  (14) 

Rather,  coupling  Is  Induced  by  the  radiation  field  associated  with  the 
dipole  operator  in  Eq.  (1), 

df^(R)  ■  <(J*1  'yl (15) 

In  the  weak-field  limit,  when  the  electric  field  strength  ISl  is  very 
small,  we  can  safely  neglect  any  higher-order  perturbations  on  the 
field-free  molecular  wavef'onctions  eind  simply  represent  the  radiative 


coupling  In  the  following  matrix  form. 


S''^'^(d)  -  -2irl  /  F"(E",R)*d(R)‘F*(E'^,R)  dR 


06) 


Obviously  this  off-diagonal,  radiatively-lnduced,  scattering  matrix 
S''^‘^(d)  has  the  same  (N“xN'^)  dimensionality  as  the  dipole  matrix  d(R). 

Because  of  the  distortions  of  the  atomic  wavefunctlons  during  the 
formation  of  the  quasi -molecule  in  Eq.  (3)  there  is,  in  general,  an 
implicit  R-dependence  Introduced  into  the  matrix  elements  evaluated  in  Eq. 
(15), 
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^uld  Eq.  (16)  can  be  evaluated  In  terms  of  two  separate  Integrals, 

srad(d)  .  s'^a^^Ca)  ♦  (18) 

The  contribution  from  the  asymptotically  vanishing  portion  b(R)»d(R)-d(») 
of  Eq.  (17)  will  always  yield  a  convergent  integral  for  the  S'*®**  (b) 
component  in  Eq.  (18).  In  fact  this  qiiantlty  is  the  source  of  the 
oscillator  strength  in  collision-induced  absorption  phenomena  (Julienne 
1982)  where  the  asymptotic  A+  hu)-+A*  transition  is  forbidden,  and  the 
asymptotic  transition  matrix  a-d(<»)  vanishes.  0<ir  concern  is  with  allowed 
transitions,  where  a*<  0,  and  we  encotmter  lindamped  oscillations  in  the 
n'jmerical  evaluation  of  the  continuum-continuum  matrix  elements  which 
makes  S'”®**  (a)  only  conditionally  convergent. 

3.  Commutator  Techniques  for  Eval»jatlng  S'*®**(a) 

Let  us  divide  the  asymptotic  elements  in  Eq.  (17)  by  the  detuning  in 
Eq.  (6)  and  obtain  the  following  dimensionless  transition  matrix  elements 
which  bear  a  direct  correlation  with  the  asymptotically  allowed  atomic 
transition  dipole, 

Ufi  ■  ^•f i^^i  (19) 

This,  in  t'irn,  is  Introduced  into  the  definition  of  the  ^'ollowing, 
Wronsklan-llke  matrix  of  continuum  wavefunctions  f'or  evaluation  at  some 


prescribed  interatomic  distance  R, 


Aa  we  shall  see,  It  is  the  asymptotic  contribution  of  V(a,<»}  to  the  integral 
srad  (a)  gq.  (^g)  and  (16)  which  causes  the  conditional  convergence  of 
the  continuum-continuum  matrix  element  in  Eq.  (1). 

Utilizing  the  following  Identity  based  on  Eq.  (10), 

f"a(T*  ♦  -  e*)r^  -  (T*  +  Y"  -  e“)P"  aF'*’  (21) 

and  noting  that  this  must  be  satisfied  for  arbitrary  values  of  the 
transition  matrix  a,  we  derive  the  following  equivalence  for  the  a 
contribution  to  the  integrand  in  Eq.  (16) 

faF*  -  F"(T"<i-or‘)F-^  +  dW(a,R)/dR  (22) 


Any  numerical  integration  over  the  first  RHS  term  is  absolutely  convergent 
because  of  the  asymptotic  decay  of  the  Interaction  matrices  V*(R)  in  Eq. 
(12).  The  complete  differential  dW/dR  merely  gives  rise  to  surface-type 
terms  which  must  be  evaluated  at  the  extremes  of  the  integration,  and,  for 
grad(j)  obtain. 


S''®‘^(a)  R*  +  X 


/  *dR  faF"^  +  ;  dR  #"(?"«-o?*)F*  +  [W(a,-)-W(a,RA)]  (23) 


We  have  chosen  to  separate  the  Integration  into  two  portions  and  only 
apply  the  equality  in  Eq.  (23)  beyond  some  arbitrarily  assigned 
interatomic  distance  R-R^*  Use  of  this  division  is  made  in  the  by  Mies  A 
Julienne  (1985)  in  the  Impact  limit  as  *  0.  For  purposes  of  this 
study,  in  the  static  wings  of  the  atomic  lineshape,  it  is  sufficient  to 
take  R^-O,  and  express  (23)  as  follows. 


S^®‘*(a)  -  -2ir  i  /  dR  P"(r<i-o?^)F*  -2iriM(a,») 


(24) 


Here  we  have  utilized  the  fact  that  both  P*  and  P'*’  are  well-behaved  radial 
wavef unctions  at  R-»0.  such  that 


M(o,0)  -  0.  (25) 

Using  the  asyoiptotic  forms  of  P^  in  Eqs.  (A1)  and  (A7)  we  obtain  the 
following  structure  for  WfjCo.R)  at  large  R, 

2Wfl(a,R)  -  I  Cafi.w^'‘^R)]S*  +  I  S*  [af.iw^"^  (R)]  (26) 

if  fi'  i'i  ff  ff  ^  f'i 

where  the  reduced  Wronskian  matrices  if^-)(R)  are  defined  in  (AS). 

Carried  out  to  infinity  we  find. 


2Wfi(a,*)  -  ^  Cafii6(e  -e  f)]S 

if  ^  f  1  i'i 


♦Is”  Caffi6(e”t-e'")]  (27) 

ff  ff»  ^  f  1  • 


As  long  as  the  photon  frequency  is  not  in  perfect  resonance  with  any  of 
the  asymptotic  atomic  transitions,  the  1t(a,»)  terms  in  Eq.  (23)  and  (24) 
will  vanish  due  to  the  delta  functions  in  Eq.  (27) .  Thus  we  can 

I 

mathematically,  and  physically,  Justify  the  neglect  of  such  non-convergent 
contributions  to  S'^^^(a)  'jnder  all  conceivable  circumstances  in  lineshape 
theory. 

If  we  use  only  single  channel  wavefiinctlons  in  Eq.  (23)  and  make  the 
choice  R^-0  then  Eq.  (24)  reduces  to  the  expression  of  Herman  and  Sando 
(1978).  If,  instead,  we  make  the  choice  that  is  large,  but  •'inite, 
beyond  some  range  of  the  V-  potentials  such  that  the  second  integral  in 
Eq.  (23)  is  negligible,  then  we  obtain  the  adiabatic  expression  of  Bienlek 
(1981).  Either  choice  of  R^^  leads  to  a  convergent  expression  for  both  the 
single  channel  and  the  multichannel  integral,  since  the  fluctuating  terms 
have  been  segregated  into  W(a,*),  identified  with  delta  functions,  and 
promptly  discarded.  We  prefer  the  choice  R^^-O,  since  this  leads  to  an 


V  s 


expression  which  Is  naturally  damped  by  the  vanishing  asymptotic 
potentials.  Thus,  the  complete  atomic  lineshape  can  be  associated  with 
the  following,  absolutely  convergent  integral, 

srad(d)  -  -2vl/*dR  #*(E",R) •Cb(R)*(f"«r’)] •r*(E‘",R)  (28) 

In  deriving  Eq.  (27)  we  have  utilized  the  fact  that  the  nuclear  angular 
momentum  quantum  numbers  If  and  In  Eq.  (a8)  must  be  equal.'  This 
is  a  consequence  of  the  asymptotic  selection  piles  for  pure  case  (e) 
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molecular  channel  states  that  are  inherent  in  Eqs.  (15)  and  (17), 

*fl  “  ®fl  ^29) 

This  constraint  does  not  apply  to  the  collision-induced  b(R)  portions  of 
the  dipole  matrix  <1(R)  and  can  in  fact  lead  to  significant  qualitative 
differences  between  allowed  and  forbidden  lineshape  behavior. 

Exactly  the  same  commutator  analysis  can  be  used  to  confirm  the 
validity  of  the  normalization  Implied  in  Eq.  (9).  If  we  replace  P"(E“) 
with  F*(E)*  and  a  with  1*  In  Eq.  (22),  we  derive  the  following 

^  d#*(E,R)*f"(E*,R)-P(E,R)*  dF*(E*,R) 

/  dR'r(E,R')  •F"^(E^,R')  -  - - -  [ -  - ] 

o  R—  2ii^g(E  -E)  dR  dR 

.  .  ,  S‘"(E)  -  S*(E*) 

-  6(E-E"^)  S*(E)  C  - ] 


(S(E-e"’) 


The  same  procedure  confirms  the  proper  normalization  of  the  outgoing 
radial  wavefunctlons  F“(E“)  in  Eq.  (9). 


nV-*, 


VV-  •.  -N  -- 
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Associated  with  each  incident  channel  is  a  radiative  number  state 

1 

ln>  which  defines  the  properties  of  the  incident  radiation  field.  A 
corresponding  number  state  |n-1>  is  associated  with  each  final  channel 
state  and  Indicates  the  absorption  of  a  single  photon  as  portrayed  in 

Eq.  (3).  In  particular,  the  electric  field  strength  in  Eq.  (15)  is 
related  to  the  n'jmber  of  photons  n  contained  in  the  source  free  volume  V 
with  the  prescribed  energy  hu  (and  other  qualities  of  interest)  as 
follows, 

l8(<ij)|2  .  2slIu(n/V)  (31) 

In  deriving  Eq.  (1)  we  have  relied  on  this  quantity  being  sufficiently 
small  to  permit  the  neglect  of  any  radiative  perturbation  of  the 

field-free  molecular  continuiom  states  y*  and  which  we  associate  with 
-  i  f 

the  bare  channels  and  respectively.  However,  for  optical 
collisions,  these  are  really  not  the  proper  scattering  channels  that 
should  be  >ised  to  define  the  absorption.  By  considering  the  proper  dressed 
atomic  states,  and  employing  the  resultant  dressed  channel  states  that  are 
generated  by  the  eisymptotlc  radiative  coupling,  we  can  gain  a  good  deal  of 
insight,  and  actually  obtain  a  more  general  perturbative  expression  for 
^rad^  To  first  order  in  E(<u)  we  will  retrieve  the  convergent 
perturbative  result  in  Eq.  (28). 


We  have  previously  developed  the  pr(^er  close-coupled  equations  for 
the  full  molecule-field  scattering  wavefunctlons  (Mies  1981) 


Oj  (n.E)  -  I^{(>*<2>ln>}  Fj,j(n,E,R)  ♦  J  (ij)”®  ln-1>}  Ffi(n,E,R)  (32) 
-  I  5j(n)Fji(n,E.R) 


Where  5j(n)  designates  an  element  on  the  N-(N'*'+M")  dimensional  row  vector 
of  molecule-field  channel  states  composed  of  products  ^^0ln>  and 
<Pj(5ln-1>  generated  from  Eq.  (7) 


t(n)  -  [••••!  ♦^®ln>,  ••••,  (>^(^n-1>,  ••••]  (33) 

The  (NxN)  matrix  of  radial  wavefunctlons  F(n,E,R)  Is  evaluated  at  the 
combined  molecule-field  energy  E 

E  -  E*  +  trflu  -  E“  +  (n-1 )  (3^) 

and  denotes  proper  Incoming  scattering  solution  vectors  originating  from 
the  following  set  of  coupled  equations. 


P(n,E,R)  -  0.  (35) 

In  the  limit,  when  n-K)  and  the  intensity-dependent  radiative  coupling 
matrix  d(n,R)  in  Eq.  (15)  vanishes,  we  retrieve  the  field-free  solution 


vectors  that  are  defined  by  Eq.  (10)  and  used  in  Eq.  (1), 


F(n,E,R) 


[F^(E^,R)  0 

0  F‘(E‘,R)  _ 


Note  that  linlike  the  molecular  interaction  matrices  V^(R)  which  are 


constructed  to  vanish  identically  as  R-»«,  the  asymptotic  radiative 
coupling  d(n,»)-*’a(n)  in  Eq.  (17)  never  vanishes  in  the  presence  of  a 


field.  Since  a  scattering  utrix  element,  such  as  (d),  is  meant  to 
describe  the  transition  probability  between  pairs  of  asymptotically 
uncoupled  states  which  can  be  independently  prepared  before  collision,  and 
asymptotically  observed  after  the  collision  event,  it  should  be  apparent 
that  our  molecule-field  channel  states  ^(n)  are  improper  scattering  states. 


Let  us  define  an  asymptotic,  orthogonal  transformation  of  oijr  channel 


states 


5D  .  p,D 


such  that  and  the  asymptotic  interaction  vanishes, 

_  rr -»*<»>  I  - H  .  (38, 

L  a(n)  T  -  e  J  L  0  T  -  w  (a)J 

The  new  set  of  diagonalized  kinetic  energies  w^(a)  (Not  to  be  confused 
with  the  reduced  Wronskian  matrices  of  the  previo«is  section  and 

Appendix  A)  obtained  from  the  secular  equation. 


j_a(n)  j  w  -  e  J 


reflect  the  dressing  of  the  asymptotic  atomic  states  associated  with  the 
transition  Aj^  +  Hflu-^Aj.*.  We  are,  of  co'irse,  concerned  with  the  weak-field 
limit  where  the  radiative  interaction  matrix  a,  and  the  corresponding 
reduced  matrix  a  in  Eq.  (19)  are  very  small.  To  second-order  we  obtain 
the  following  estimates  for  the  field-shifted  asymptotic  kinetic  energies 
of  the  incident,  and  final  molecule-field  channel  states 


w'*‘(a)  -  e*  -  (30)44  -  e*  -  J]  a  /(e~-e*) 
i  i  i  f  fi  f  i 


(MOa) 


w  (a)  -  e  “  (aa)ff  -  e  -  /  a  /(e  -e  ) 


fl  f  1 


(40b) 


If  we  recognize  that  a^j  13  related  to  the  Rabl  frequency  Qf^  for  the 
asymptotic  atomic  transition  fif^saafi^/h  and  note  that  the  atomic  detuning 
frequency  (6u)f^  Is  related  to  the  asymptotic  kinetic  energy  difference  In 
Eq.  (6),  (6u}f^  «  Afi^/h,  we  find  that  the  shifts  In  Eq.  (41)  are  Identical 
to  the  usual  expressions  for  dressed  atomic  states  in  we2bc-fleld  limit 
(Mies  and  Ben  Aryeh  1981,  Mies  1983). 


Vi  [1 


AfiAfu 


-  ^1  Cl  •*fl] 


(41a) 


(4lb) 


Let  us  define  a  field-shifted  reduced  matrix  x^ comparable  to  cif^  in 


Eq.  (19), 


w  -  w 
f  1 


«fl^1"''fi) 


To  second-order  in  z  the  orthogonal  matrix  may  be  expressed  as 
follows. 


U».  P* : 1 . .  - 1 

[  z  I  1  -  z  z/2  J 


Transforming  the  expansion  In  Eq.  (32)  to  the  dressed-channel  basis 
a"’  -  CF(n,S,R)  -  (EU°)  (0°F(n,E,R))  -  C°F°(n,E,R)  (44) 

we  obtain  the  following  set  of  coupled  equations  for  the  radial 
wavef'inctions  F® 


rj*  *  M*  -  w'*’(a)  ♦  Y*  b  -  (V"’z-z¥  )  1  „ 

-  F°(n,E, 

b  +  (ITz-zr)  r  +  V"  -w~(a)  ♦  I" 


R)  -  0  (45) 


r:  rrr?: 

>  *  >  V  v" 

J'"'. 


•v 

‘V 

'V  'i-v: 
'""is 

,  .  .'.sV 


where  the  T^(R)  matrices  vanish  as  R-**  and  only  contain  second-order  terms 


I*(R)  -  i  T"x  >  (Sx+fe)  (46a) 

I"(R)  -  *  Y+i  -(xx¥"+V"xx)/2  -  (x6+bi)  (46b) 

In  the  first-order  limit  we  can  see  that  the  perturbative  solution  to  Eq. 

(45)  yields  the  same  diagonal  blocks  of  field-free  solutions  as  Eq.  (36), 
with  an  off-diagonal  radiative  coupling  b-*-(T“«raY*’)  matrix  which  yields 
the  approximation  (Taylor  1972,  Mies  1981) 

S>'acl(d)  :  -2iri  <F"(E")*lb  T"<r-a?*l  P*(E‘*’»  (47) 

In  complete  agreement  with  our  commutator  analysis  in  Eq.  (28). 

Obviously  by  referring  our  radiative  scattering  matrix  elements 
to  the  proper  dressed  channels  and  rather  than  Cf  and  we  have 
performed  the  equivalent  of  ignoring  the  indeterminant  surface  integral 
M(o,®)  in  Eq.  (24).  In  deriving  Eq.  (47)  we  have  neglected  the  shift 
parameter  in  Eqs.  (41)  and  (42)  and  taken  z  -  a.  However,  the 
close-coupling  analysis  is  not  committed  to  small  perturbation,  and  we  can 
progress  to  second-order  by  solving  Eq.  (45),  or  to  arbitrary  order  in  the 
radiative  coupling  if  we  solve  the  coupled  equations  (35)  directly, 
subject  to  the  dressed  state  boundary  conditions  implied  by  Eq.  (44).  For 
example,  the  dressing  of  an  isolated  atomic  transition  always  insures  that 
^fi  £  *fi  £  ®fl*  effect  introduces  a  saturation  of  the 

pressure-broadened  half-width  in  intense  laser  fields  (Mies  and  Julienne 
1985). 
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5.  Dlscusalon 


The  numerical  solution  of  the  close-coupled  equations  (45)  between 
dressed  levels  never  causes  any  difficulties  associated  with  conditional 
convergence.  There  Is  an  absolute  convergence  of  the  radiative 
scattering  matrix  element  obtained  from  imposing  the  following 
asymptotic  structure  on  F*^(n,E,R) 

F°  (n,E,R)  -  (i/2)C6fi  fJ"^w"’(a),R)  -  gi^\w‘(a).R)  (48) 

f  1  R+«  11  f  f  f  1 

whe.'e  the  appropriate  Incoming  f|”^  and  outgoing  g^'*'^  Bessel  f'jnctions 
are  defined  in  Appendix  A.  However,  care  must  be  taken  to  use  the 
correctly  shifted  asymptotic  kinetic  energies  defined  by  the  secular  Eq. 
(39)  rather  than  the  field-free  values  and  .  Essentially  what  we 
are  defining  In  Eq.  (48)  are  the  particular  set  of  degenerate  continuum 
wavef'inctiona  needed  to  describe  the  collision-induced  absorption  between 
the  stationary  field-dressed  states  of  the  asymptotic  atom  A  in  the 
radiative  field  associated  with  the  Rabi  frequency  in  Eq.  (42). 

This  close-coupled  procedure  avoids  any  commitment  to  weak  fields  and 
perturbative  estimates  of  .  However,  to  first  order  in  the 

asymptotic  radiative  coupling  af^  the  dressed-state  eigenvalues  in  Eq. 

(41)  are  unshifted,  and  (d)  is  given  by  the  first-order  distorted 

wave  approximation  in  Eq  (47).  By  referencing  O'lr  close-coupled 
wave fijnct ions  to  the  dressed  channels  and  rather  than  and 
in  Eq.  (58)  we  are  essentially  performing  the  equivalent  of  ignoring  the 
indeterminant  surface  Integfral  W^^(a,*)  in  Eq.  (24). 


We  have  performed  multichannel  close-coupled  calculations  of  the 
Sr(^S)  Sr(^P)  transition  broadened  by  collision  with  Ar  (Julienne,  1983, 
1985).  One  component  of  the  absorption  rigorously  reduces  to  a  simple 
radiative  transition  between  a  single  channel,  elastically  scattered, 
ground  state  continuum  T*  (E*)  -  with  a  field-free  radial 

function  described  by  scattering  from  the  isolated  molecular  potential 
V*^(R)  with  a  resultant  phase  shift  n^CEj) 


11 


(E*,R) 


sin  (kjR— itli/2+ni) 


(49) 


and  a  single  channel,  elastically  scattered,  excited  state  contlnuiim 
^^(E")  ■  *  phase  shift  nf(ep  determined  by  scattering 

off  the  molecular  potential  V~^(R) 


F"^  (E“,R)*  -  e  ain(kfR  -ir4^/2  ♦  nf)^  (50) 

In  the  spectroscopic  terminology  one  would  associate  with  Eq.  (3),  we  are 
considering  the  Q-branch  continmam-continu'm  absorption  of  the  ^  E(X)-*-^  n(A) 
transition  of  the  (SrAr)  quasi molecule. 

A  comparison  of  the  exact  two  state  close-coupling  matrix  elements 
obtained  from  Eq.  (48)  and  the  numerically  evaluated  matrix  elements  in 
Eq.  (47)  confirms  the  weak -fie Id  result 


rad  -  - 

S  (a)  -  -2iri  a^,,  <F  (E“) 

fi  a-K)  ff 


(51) 


This  is  demonstrated  in  Table  1  for  a  range  of  detunings.  Herman  and 
Sando  (1978)  have  already  derived  and  made  use  of  this  limited  two-state 
version  of  Eq.  (47)  for  calculating  the  wings  of  resonant  Li  and  Na  lines 


This  expression  allows  for  a  convergent,  fully  quantal  evaluation  of 
optical  collision  cross-sections  without  recourse  to  the  usual,  although 
generally  quite  satisfactory  (Mies  1968,  Jablonskl  1945)  expressions 
based  on  the  JWKB  and  stationary  phase  a|>proxlfflatlons  (see  Allard  and 
Kielkopf  1982).  Although  we  have  made  a  more  explicit  analysis  of  the 
Indeterminant  character  of  the  matrix  element  (a),  which  requires  the 
isolation  of  the  asymptotic  surface  term  in  Eq.  (24),  we  conclude  that 
this  can  legitimately  be  ignored  from  f«irther  consideration.  Thus  we  have 
used  a  more  rigorous  and  general  procedure  to  Justify  the  specialized 
Herman  and  Sando  result  in  Eq.  (51). 

Our  more  general  expression  in  Eq.  (47)  does  Include  multichannel 
couplings  and  allows  for  the  possibility  of  inelastic,  and  non-adiabatic 
contributions  to  atomic  llneshapes.  Again,  a  comparable  version  of  this 
multichannel,  first-order  expression  has  appeared  in  the  llneshape 
literature,  although  it  has  been  specialized  to  the  Stark-broadening  of 
the  hydrogen  atom  (Feautrler  et  al  1976).  A  more  general  multichannel 
expression  for  Stark-broadening  can  be  extracted  from  Seaton's  (1981) 
analysis  of  free-free  transitions  in  electron-atom  collisions,  although  we 
have  made  a  more  explicit  outline  of  the  proper  bo<indary  conditions,  and  a 
less  restrictive  application  of  the  theory. 

We  would  claim  that  our  major  contribution  is  to  recognize  that  the 
convergent  first-order  matrix  element  in  Eq.  (47)  corresponds  to  the 
scattering  among  asymptotic  dressed  states  of  the  absorbing  atom.  Since 
the  dressed  field-atom  states  in  Eq.  (37)  are  defined  as  stationary 


eigenfunctions  of  the  complete  fflolecule'*f ield  Hamiltonian  at  infinite 
interatcmic  separation  we  have  removed  the  troublesome  coupling  in  the 
asympototic  regions  of  the  molecular  continuum  states. 

part 

This  finding  means  that  we  can  now  confidently  equate  in  Eq.  (1) 
with  the  numerically  derived  elements  obtained  from  solution  of  the 
close-coupled  equations  (35)  transformed  to  the  appropriate  boundary 
conditions  in  Eq.  (48) .  In  fact,  since  stable  close-coupling  codes  are 
required  to  generate  the  multichannel  wavef'jnctions  In  Eqs.  (8)  and  (10), 
It  is  often  a  simple  matter  to  merely  expand  the  set  of  coupled  equations, 
as  in  Eq.  (35)  or  (45),  to  directly  Include  the  radiative  coupling  between 
channel  states  as  required  by  the  boundary  conditions  in  (48)  and 
completely  avoid  explicit  numerical  quadrature  of  the  multichannel  matrix 
elements  in  Eq.  (47).  Obviously,  this  convenience  Is  accomplished  by 
sacrificing  computer  time  to  solve  for  larger  sets  of  coupled  radical 
equations. 

When  performing  the  close-coupled  calculations  aoxjng  dressed  levels 
r*3cl  I  I 

to  extract  S^^'^(d)  we  must  choose  a  representative  field  strength  1b((ii)I 

part 

in  Eq.  (15)  sufficiently  small  to  insure  that  (d)  in  Eq.  (48) 
remains  linearly  proportional  to  the  radiative  coupling  strength  d(R)  in 
Eq.  (35)  and  yet  sufficiently  large  to  ins'ire  the  numerical  accuracy  the 
resultant  scattering  matrix  elements.  As  long  as  the  magnitude  of  (d) 
remains  below  about  10~3  it  seems  to  remain  quite  linear  in  d,  and  there 
is  a  one-to-one  correspondence  between  Eq.  (28)  and  the  values  extracted 
from  (48).  This  is  substantiated  by  the  results  in  Table  2.  Of  coiirse 
tnie  deviations  from  linearity  must  ultimately  prevail  as  the  actual  field 


strength  in  Eq.  (31)  Increases,  reflecting  the  onset  of  saturation  of  the 
contlnuiffli-continu’im  transition  in  Eq.  (3).  This  is  most  prominent  near 
line-center,  and  is  given  greater  attention  in  the  forthcoming  paper  on 
the  impact  limit  (Mies  and  Julienne  1985). 
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Appendix  A;  Aaymptotlc  Properties  of  Multichannel  Contlnmam  Wavef'inctlons 

In  contrast  to  conventional  scattering  notation  we  have  <jsed  *  and  - 

to  label  two  separate,  mutually  orthogonal,  sets  of  Initial  and  final 

channel  states  each  with  Its  own  matrix  of  close-coupled  radial 

solution  vectors  P*  and  P"  respectively.  These  field-free  scattering 

wavefunctlons  are  expected  to  adequately,  and  independently,  describe  any 

inelastic  couplings  among  the  set  of  initial,  and  the  set  of  final  atomic 

states  Involved  in  the  absorption.  Both  F*  and  P~  are  assigned  the  usual 

incoming  boundary  conditions  appropriate  to  defining  the  inelastic 

cross-sections  Induced  by  collisions  with  the  perturbing  atom.  Once 

radiative  coupling  is  introduced  between  the  two  seta  of  channels  we 

associate  with  the  subset  of  incoming  channel  states  defined  by  the 

absorption,  while  becomes  identified  with  the  subset  of  excited 

outgoing  channel  states.  The  conventional  boundary  conditions  in  F'^  are 

♦ 

adequate  to  define  the  incoming  contlnu'im  state  required  in  Eq.  (1). 

The  specific  normalization  of  P*  and  P“  that  we  introduce  in  this  appendix 
is  carefully  chosen  to  insure  that  (P“)*  also  defines  the  proper  outgoing 
scattering  state  in  Eq.  (8b)  and  (1). 


The  radial  wavefunctlons  in  Eq.  (8a)  have  the  following  asymptotic 
structure 

P'’(E'",R)  -  i/2  Cf^~^€'^,R)  -  f^'’^(e'",R)  S'’(E"’)]  (AT) 

where  f^^^(e'*',R)  are  diagonal  (N*  x  N'*’)  matrices  of  energy-normalized 
incoming  (-)  and  outgoing  (+)  spherical  waves  [see  Roman  (1965)  for  the 
distinction  between  Incoming  f'*'  and  outgoing  (T^)*  scattering  states  and 
these  spherical  Bessel  ^unctions], 


f  (e^.R) 
i  1 


n-*» 


exp[±i(kjR-4j/2)] 


(A2) 


Each  incoming  channel  state  has  a  well-defined  asymptotic  kinetic 
energy  e*  which  is  prescribed  by  the  total  energy  S'**,  which,  in  turn, 
determines  the  magnitude  of  the  asymptotic  radial  wavevector  k^ 

kJ(E*)  -  2uj^BeJ/A2  -  (A3) 

The  asymptotic  constraint  imposed  on  the  interaction  matrix  V'*'(R)  in  Eq. 
(12)  insures  that  each  of  the  resultant  case  (e)  channel  states 
correlate  with  a  specific  nuclear  angular  momentum  quantum  number,  or 
partial  wave  These  running  waves  only  occur  for  open  channels  where 
k^(E'*')  is  real  and  positive,  such  that 
(±)'  (±)  it) 

f  S  df  /dR  -  ±ik.  f  (A4) 

1  i  R*«  ^  i 


The  particular  normalization  in  Eq.  (A2)  insures  the  following  asymptotic 
property 


\cfr  fj-  fj-  <T'i 


(±)  (T)’  (±)'  (T) 


(A5) 


(ei  -  ej) 


This  is  rationalized  in  delta  function  theory  by  appealing  to  the 
subsequent  effect  of  averaging  of  Eq.  (A6)  over  small  increments  of 


energy  about  (e^  -  tj ) .  This  procedure  is  physically  justified  because 


scattering  eigenstates,  which  are  not  normalizable,  are  ultimately  used  to 


expand  states  for  which  the  square  integrable  normalization  does  exist. 
(See  Messiah  (1958)  where  this  procedure  is  discussed  in  detail  and  its 


equivalence  to  other  methods  of  generalizing  orthonormality  relations  to 
the  continuum  is  demonstrated) .  The  saiM  limiting  procedure  is  used  to 
show  that  we  may  properly  write. 


.a)  0 


(A6) 


A  structure  comparable  to  ?*(£*, R)  la  defined  for  incoming  radial  solution 


vectors  F"(E“,R)  associated  with  the  set  of  outgoing  channels  in  Eq.  (10), 

f 


r(E',R)  -  i/2Cg(")(c‘,R)  -  g^-^>(e-,R)  S"(E")] 

R+« 


(A7) 


In  place  of  f^*^(e‘*',R)  we  introduce  diagonal  (N“  x  N~)  matrices  of 
spherical  Bessel  functions  g^^^(e“,R)  as  defined  by  the  asymptotic  kinetic 
energy  e“,  wavevector  kf  and  partial  wave  quantum  number  If  associated 
with  channel  state  Actually  these  nonning  waves  exhibit  the  same 
structure  and  properties  given  for  f^*^  in  Eqs.  (A2)-(A6).  In  fact  the 
cross-terms  between  these  sets  will  also  conform  to  Eqs.  (A5)  and  (A6) 
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However,  we  have  introduced  the  distinct  notation  g  to  emphasis  that 
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these  matrices  as  compared  to  f'*'  are  of  different  dimension,  are 


associated  with  a  separate  sub-set  of  channel  states,  and  are  evaluated  at 
a  different  total  energy  E“  as  prescribed  by  Eq.  (5). 
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Solving  the  two  sets  of  coupled  Eqs.  (10)  separately  at  total 
energies  E'*’  and  E"  we  obtain  an  independent  pair  of  symoetric,  unitary 
scattering  matrices  S'*’(E'*’)  and  S~(E~)  which  separately  embody  all  the 
inelastic,  non-adlabatlc  couplings  among  the  incoming  and  outgoing 
channel  states  respectively, 

S±(S±)*  -  SHS^)*  -  (A10) 

The  boundary  conditions  Imposed  on  Eqs.  (A1 }  and  (A7)  have  been  devised  to 
Insure  that 

P^(E±,R)  -  P±(E±,R)*  S±(E±)  (All) 

which  la  necessary  to  satisfy  Eq.  (13).  Further,  given  that  both  F*  and 
F“  satisfy  proper  boiindary  conditions  for  the  usual,  physically 
significant,  incoming  states  of  scattering  theory,  the  result  in  Eq.  (All) 
confirms  that  the  complex  conjugate  function  P“(E“,R)*  we  have  used  in  Eq. 
(8b)  does,  in  fact,  represent  the  proper  outgoing  state  that  we  require  in 
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Table  I.  Conparlson  of  cloae-coupled  and  perturbative  radiative 

scattering  matrix  elements  evaluated  for  radiative  coupling 
strength  a-1 .674x10~‘*  cm."^  corresponding  to  a  radiative  flux 
F„  -  10“^  watts/cm2  in  Sr  ♦  Ar  system. 


1  radi 
'Si2  I 

.  '"ad  V, 

arg(Si2  )° 

1  radi 

>Si2  1 

,  rad 
arg(Si2  ) 

Is^fl 

,  rad 
arg(Si2 

-50.0 

1.1979(-5)° 

1.1979(-5)<^ 

2.1979 

2.1976 

1.4086(-5) 

1.4086(-5) 

-0.1151 

-0.1153 

2.42l8(-5) 

2.4217(-5) 

1.8389 

1.8371 

-20.0 

4.3643(-6) 

4.35aa(-6) 

2.7461 

2.7457 

1.3909(-5) 
1.3911 (-5) 

0.6519 

0.6518 

3.3967(-6) 

3.3959(-6) 

-0.8593 

-0.8611 

-  1.0 

3.4027(-4) 

3.4027(-4) 

1 .1641 
1.1637 

2.8557(-4) 

2.8558(-4) 

-0.8012 

-0.8013 

3.3916(-4) 
3.391 4(-4) 

-0.6293 

-0.6311 

+  0.1 

3.2451(-3) 

3.2451(-3) 

-2.0678 

-2.0681 

2.5898(-3) 

2.5399(-3) 

2.2575 

2.2573 

3.1094(-3) 

3.i097{-3) 

-0.6171 

-0.6190 

♦  1.0 

3.1020(-4) 

3.1019(-4) 

-2.1416 

-2.1418 

2.3674(-4) 

2.3678(-4) 

2.1896 

2.1895 

2.8671 (-4) 
2.8675(-4) 

-0.6074 

-0.6092 

♦20.0 

1.2597(-6) 

1 .2576(-6) 

-0.5339 

-0.5343 

5.3994(-6) 

5.3982(-6) 

-2.3647 

-2.3649 

2.9928(-6) 

2.9930(-6) 

-0.4199 

-0.4218 

♦50.0 

1.29a6(-6) 

1.29a5(-6) 

-2.8747 

-2.8750 

5. 9208 (-7) 
5.9233(-7) 

1.7570 

1.7568 

3.7391 (-8) 

3.7342(-8) 

-0.1986 

-0.2006 

♦100.0 

1 .2478(-7) 

1 .2459(-7) 

2.8492 

2.8489 

3.5848(-8) 

3.5828(-8) 

1.4673 

1.4671 

6.5553(-9) 

6.5774(-9) 

-0.0158 

-0.0177 

(a)  J  Is  the  total  angular  momentum  of  the  Initial  and  final  molecular  continuum 
states  in  the  quasi molecular  Sr+Ar  Q-bramch  transition. 

pari 

(b)  Argument  of  the  complex  S^^  tnatrix  element  In  radians. 

(c)  Upper  element  In  each  pair  Is  obtained  from  solution  of  coupled  equations. 


(d)  Lower  element  in  each  pair  of  entries  is  evaluated  by  numerical  quadrature  of 
convergent  continuum-continuum  matrix  element  in  Eq.  (51). 


Table  II.  Variation  of  close-coupled  radiative  scattering  matrix  element  with 
coupling  strength. 


2  -  1 .0  cm 


arg(Si2  )  SioA 


2  ■  “50.0  cm 


arg(Si2  )  Si9A 


1.0(-6) 

1.0(-4) 

1.0(-2) 

1.0(+0) 

1.0(+2) 

1.0(+4) 

1.0(+6) 

1.0(+8) 

1.0(+9) 

1.0(+10) 


(a)  radiation  flux  in  units  of  watts/cm^.  Using  predicted  transition  dipoles  for 
Sr+Ar  this  flux  yields  the  listed  radiative  coupling  coefficient 


1.67(-6) 

3.4027(-6) 

1.1641 

0.32352 
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2.1979 
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0.39441  ^ 
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0.00000 
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Growth  and  Decomposition  of  Graphitic  Structures 

Principal  Investigators:  R.  L.  Brown  and  S.  E.  Stein 

Introduction 

Graphitic  materials  are  the  most  commonly  discussed  laser-armor 
protective  materials  for  spacecraft.  Upon  laser  heating,  the  molecular 
structures  in  these  materials  decompose  and  reform,  probably  in  both  the 
vapor  and  condensed  phases.  At  present  there  exists  no  molecular  under¬ 
standing  of  these  chemical  processes.  This  progreun  Intends  to  develop 
this  understanding  in  order  to  assist  both  the  rational  analysis  of 
laser-induced  ablation  processes  as  well  as  the  development  of  improved 
protective  materials. 

The  nature  and  reactivity  of  the  chemical  structures  formed  in  these 
processes  can  control  the  effectiveness  of  laser  protective  materials  in 
several  ways.  (1)  The  formation  rate  and  nature  of  the  volatiles  produced 
may  depend  on  both  the  nature  and  concentrations  of  precursor  sites  in  the 
solid  as  well  as  on  the  reaction  mechanism.  (2)  The  chemical  structure 
and  reactivity  of  the  materials  remaining  after  ablation  by  the  initial 
laser  pulse  will  be  determined  largely  by  chemical  factors.  (3)  The 
condensation  of  carbon  vapors  to  form  refractory  particles  may  not  be 
controlled  by  diffusion,  but  by  chemical  reaction,  since  reformed  struc¬ 
tures  must  form  strong  chemical  links  to  resist  re-evaporation. 

The  development  of  realistic  working  models  for  the  above  reactions 
requires  a  substantial  increase  in  our  understanding  of  the  structure  and 
reactivity  of  graphitic  species.  Our  combined  theoretical/experimental 
program  focuses  on  the  fundamental  chemical  properties  of  very  large 


polyaromatlc  structures.  Graphite  and  nearly  all  carbonized  organic 
materials  are  composed  primarily  of  these  structures. 

Early  In  this  program  we  discovered  that  the  atomic  orientation  at 
the  edges  of  these  structures  and  their  overall  size  were  critically 
important  factors  in  determining  their  reactivity.  Our  theoretical  work 
has  therefore  been  focused  on  separating  and  examining  Individually  these 
two  variables.  Our  experimental  program  Is  devoted  to  examining  well 
defined  carbon  deposition  reactions  which  can  be  connected  with  our 
theoretical  work. 


During  this  period  we  used  the  Pople  self -consistent-field  (SCF) 
molecular  orbital  theory  to  calculate  orbital  energy  levels  and  pi- 
electron  distributions  of  large  hexagonal  polyaromatlc  hydrocarbons.  As 
in  our  previous  work,  four  series  of  molecules,  each  with  a  different  edge 
structure  was  considered.  For  molecules  having  an  anthracene-like  edge 
structure,  SCF  calculations  predict  that  the  pi-electrons  in  the  highest 
occupied  orbital  (HOMO)  will  be  concentrated  at  the  perimeter  of  the 
molecule.  These  high  level  electrons  are  the  ones  most  likely  to  become 
Involved  in  reactions.  This  concentration  of  electron  probability  at  the 
edges  of  the  molecule  is  not  limited  to  the  highest  level.  It  is  evident 
in  orbitals  whose  energies  lie  as  much  as  0.8  eV  below  the  highest  level. 
In  contrast  to  these  molecules  are  those  with  phenanthrene-like  edges. 

For  them,  the  electrons  In  these  high  orbitals  are  uniformly  distributed 
over  the  whole  molecule.  Because  of  the  large  number  of  orbitals  in  these 
big  molecules  it  was  not  feasible  to  examine  the  pi-electron  distribution 


In  each  orbital.  Instead  we  calculated  the  average  position  of  an 
electron  In  a  p^ticular  orbital.  From  this,  those  orbitals  having 
unusual  electron  distributions  could  easily  be  identified. 

We  next  applied  the  simpler  Huckel  molecular  orbital  (HMO)  theory  In 
a  systematic  way  to  the  oiolecules  In  our  different  homologous  series. 

Like  SCF  calculations,  HMO  calculations  yield  orbital  energies  and  wave 
functions  for  pl-electrons .  From  these,  a  number  of  molecular  properties 
may  be  calculated.  So  far  we  have  examined  energy  level  densities,  bond 
orders,  electron  distributions,  free  valence  (a  reactivity  Index)  and 
resonance  energies.  Over  thirty  years  ago,  Coulson  and  Taylor  [1]  applied 
HMO  theory  to  graphite  considered  as  an  infinite  planar  polybenzenold 
molecule.  They  calculated  the  pl-electron  energy  per  carbon,  the  density 
of  states  In  the  pi-electron  bond,  and  pi-bond  order.  We  have  been  able 
to  compare  our  calculations  of  the  properties  for  large  finite-sized  mole¬ 
cules  with  their  limiting  values  for  graphite. 

The  density  of  pi-state  levels  as  a  function  of  energy  is  shown  in 
Figure  1  for  a  very  large  molecule  having  281 4  carbons  and  a  phenanthrene 
type  edge.  The  smooth  curve  shows  the  graphite  limit  calculated  by 
Coulson  cuid  Taylor.  The  density  fluctuations  about  the  graphite  limit 
observed  at  the  lower  energies  are  expected  to  damp  out  in  larger 
!ix)lecules. 

For  a  smaller  molecule  (1014  carbons)  with  this  same  edge  type. 

Figure  2  shows  how  the  pi-bond  order  varies  as  one  moves  from  the  center 
to  the  perimeter.  Only  the  bond  orders  for  bonds  near  the  edges  of  the 
molecule  differ  appreciably  from  the  graphite  limit  shown  by  the  dotted 


The  orders  of  the  bonds  in  the  centers  of  our  molecules  are  given  in 
Figure  3  as  a  function  of  molecular  size  (Nc  -  number  of  carbons) .  They 
clearly  extrapolate  to  the  graphite  limit  in  each  of  the  different  series. 

Hess  and  Schaad  [2]  have  proposed  a  resonance  energy  calculated  as 
the  difference  between  the  Huckel  pi-electron  energy  and  that  of  a 
particular  reference  structure  which  they  define.  Figure  4  shows  the 
resonance  energies  per  carbon  for  each  of  these  series  as  a  function  of 
molecular  size.  They  all  extrapolate  to  the  graphite  limit  although  the 
series  containing  anthracene  type  edges  (series)  experiences  a  deep 
minimum  in  resonance  energy  at  about  400  carbons. 

In  graphite,  the  energy  of  the  highest  occupied  molecular  orbital  is 
at  the  energy  zero.  For  all  of  the  series,  the  energies  of  the  HOMO'S 
extrapolate  to  zero.  However,  as  shown  in  Figure  5,  the  rate  of  approach 
to  this  limit  varies  widely.  Molecules  with  phananthrene  type  edges 
(series  3)  have  HOMO  levels  which  are  particularly  Slow  to  approach  zero. 
Since  the  pi-electrons  in  these  top  levels  are  the  ones  most  likely  to  be 
involved  in  chemical  reactions,  this  suggests  that  molecules  having  edges 
of  the  phananthrene  type  would  relatively  non-reacttve  until  they  reached 
a  very  large  size.  A  similar  prediction  regarding  reactivity  is  shown  in 
Figure  6  where  the  maximum  free  valence  is  plotted  as  a  function  of 
molecular  size. 

As  in  the  SCF  results,  these  variations  between  the  HOMO'S  of  the 
different  series  are  also  quite  evident  in  the  pi-electron  distributions. 
In  the  more  reactive  anthracene  edge  series,  electrons  in  this  top  level 
are  concentrated  around  the  molecular  perimeter,  while  in  the  pharanthrene 
edge  aeries  these  electrons  are  distributed  uniformly  over  the  whole 
molecule.  This  is  shown  in  Figures  7  and  3  for  molecules  of  comparable 


size  having  these  edge  types.  In  the  anthracene  edge  type,  this  concen¬ 
tration  of  electron  density  at  the  edge  is  not  restricted  to  the  top  level 
only  but  extends  downward  some  seven  or  so  levels. 

Of  all  the  theoretical  methods  that  we  have  examined,  HMO  provides 
the  most  satisfactory  all-around  approach  to  the  calculation  of  a  number 
of  properties  of  large  polyaromatic  molecules.  It  has  been  quite  success¬ 
ful  in  predicting  electronic  properties  of  small  benzenold  molecules. 

With  an  appropriate  choice  of  adjustable  parameters,  Hess  and  Schaad  have 
shown  that  it  can  yield  very  accurate  predictions  of  aromaticity  [2]  and 
heats  of  formation  [3]  for  a  variety  of  compounds.  These  workers  did  not 
include  graphite  or  steric  interactions  in  their  parameterization  scheme. 
Preliminary  work  by  us  has  indicated  that  their  inclusion  will  improve  the 
accuracy  of  the  HMO  method.  Because  HMO  predictions  extrapolate  to  the 
correct  graphite  limits,  we  feel  that  this  theory  will  also  accurately 
predict  properties  of  intermediate  sized  polybenzenoid  molecules. 
Consequently,  it  promises  to  be  a  significant  aid  in  devising  models  for 
reactions  of  carbonaceous  substances. 
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Experimental 

We  have  examined  the  kinetics  of  phenyl  radical  (  )  deposition 

on  pyrolytic  graphite  surfaces.  These  radicals  were  generated  by  the 
thermal  decomposition  of  nitrosobenzene,  CgHg-NO,  in  a  very  low  pressure 
flow  reactor.  The  pyrolytic  graphite  was  formed  by  the  deposition 
reaction  itself  on  the  fused  silica  reactor  surface.  The  key  result  of 
these  experiments  is  that  the  probability  of  deposition  per  molecule 
depends  on  the  square  root  of  the  concentration  of  radicals.  This  is 
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interpreted  as  resulting  from  an  increase  in  surface  active  sites  with 
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Increasing  deposition  rate.  Specifically,  assuming  that  the  probability 
for  deposition  per  collision  with  a  reactive  site  is  constant,  the  density 
of  reactive  sites  is  proportional  to  the  square  root  of  phenyl  radical 
concentration. 

The  probability  of  deposition  per  collision  was  of  the  order  of  1/100 
to  1/1000  per  collision  in  our  experiments.  This  suggests  that  the 
exposed  pyrolytic  graphite  surface  is  quite  reactive  and  probably  disor¬ 
dered;  this  is  in  sharp  contrast  to  the  known  crystallographic  regularity 
of  bulk  pyrolytic  graphite. 

Several  potentially  reactive  hydrocarbon  molecules  have  been  examined 
for  their  surface  reactivity  in  our  reactor.  Unfortunately,  we  have  not 
found  one  that  underwent  significant  deposition  over  our  currently 
accessible  temperature  range  (<1300K).  It  is  now  clear  that  to  examine 
deposition  kinetics  in  detail,  we  will  have  to  modify  our  apparatus  to 
enable  operation  at  significantly  higher  temperatures  (i.e.,  above  the 
melting  point  of  fused  silica).  Modifications  are  in  progress. 
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Figure  Captions 

Figure  1 .  Density  of  HMO  orbital  energy  levels  as  a  function  of  energy  for 
molecule  with  28l4  carbons  and  a  phenanthrene  edge  type.  The 


smooth  curve  shows  the  graphite  limit. 


Figure  2.  Variation  of  pi-bond  order  throughout  a  molecule  with  1014 

carbons  and  phenanthrene  type  edge.  Dotted  line  shows  graphite 
value . 


Figure  3.  3ond  orders  of  bonds  in  centers  of  molecules  in  different 

homologous  series  as  a  function  of  molecular  size.  Dotted  line 
shows  graphite  limit. 


Figure  4.  Resonance  energies  per  carbon  as  a  function  of  molecular  size. 


Dotted  line  shows  graphite  limit 
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Figure  5.  Energy  of  highest  occupied  noleculcu*  orbital  as  a  function  of 
axjlecular  size. 


Figure  6.  Maximum  free  valence  as  a  function  of  molecular  size. 
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Figure  7.  Electron  probability  distribution  .of  highest  occupied  molecular 
orbital  in  a  molecule  with  anthracene  edge  type. 

Figure  8.  Electron  probability  distribution  of  highest  occupied  molecular 
obital  in  a  molecule  with  phenanthrene  edge  type. 
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Flgur*  2.  Variation  of  pl-bond  order  throughout  a  aolacula  with  1014 

earbona  and  phananthrana  type  adga.  Dotted  line  shows  graphite 
value. 


Figura  3.  Bond  ordors  of  bonds  in  centers  of  oolecules  In  different 


hOBologeus  series  as  a  function  of  solecular  size.  Dotted  line 
snows  grapttlte  Halt. 
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BUILD-UP  AND  I  ERADIATION  OF  OBSCURIEK3  CLOUDS  UNDER  NEAR-VACUUM  COEDITIONS: 

APPLICATION  TO  SPACECRAFT  SURVIVABILITY 

Principal  Investigators:  Hwmas  J.  Buckley 

Sharon  G.  Lias 

Annual  E^eport 
FY  1985 

This  project  concentrates  on  the  elucidation  of  the  rates  and  mechanisms 
of  ionic  polymerization  mechanisms  in  highly  unsaturated  organic  conpounds, 
the  structures  and  thermochemistry  of  relevant  ions,  the  photodissociation  of 
the  polymeric  ions,  and  the  repair  mechanisms  which  regenerated  the  polymeric 
ions  after  photodissociation.  Positive  ions  generated  in  a  pulsed  ion 
cyclotron  resonance  spectroneter  undergo  rapid  polymeriza-ion  with  surrounding 
neutral  molecules  (v^ic±  can  be  Identified  with  molecules  vaporized  from  a 
graphitic  material  at  the  surface  of  a  spacecraft) .  In  the  kinetic  mass 
spectrcmeters  available  for  this  project  it  is  possible  to  reproduce  some  of 
the  conditions  whicii  would  prevail  near  the  surface  of  a  spacecraft  both  in 
terms  of  gas  pressure  and  composition,  and  in  terms  of  impinging  natural 
electromagnetic  and  laser  radiation. 

In  investigations  of  the  generation  and  dispersal  of  smokes  and  obscuring 
materials  in  space,  the  kinetics  of  evaporation  of  refractory  materials  under 
laser  irradiation,  and  condensation  processes  leading  to  the  formation  of 
microscopic  particles,  graphitic  materials  and  polyary lacetylenic  polymers  are 
of  particular  interest.  Therefore,  the  work  completed  during  FY85  focussei 
pci  aril/  on  the  elucidation  of  the  gas  phase  polymerization  and  reaction 
mechanisms  following  laser-irradiation  induced  breakdown  in  aromatic 
hydrocarbons,  particularly  those  with  acetylenic  and  olefinic  side  chains 
such  as  phenylacetylene  and  phenylethylene  (styrene) .  In  addition,  the 
kinetics  and  mechanisms  of  the  reactions  of  phenyl  ions  were  examined,  and 


structures  of  CgHg'*’  ions  generated  from  a  number  of  organic  molecules  have 
been  elucidated. 


Pheny lacety lene 

The  ion-rnolecule  chemistry  in  the  pheny lacety lene  system  has  never  before 
been  investigated.  Our  experiments  showed  that  the  predominant  reactions  of 
the  parent  ion  with  the  pheny lacety lene  molecule  are: 


^Rn 

an'^/molecule-sec 

rtlO 


CgHg^  +  CgHgCaCH  ->  Products 

■*>  ‘^16^12^ 

V  LI  ▼ 


9.1  +0.5 


^16^11+  ^ 

^16^10+  ”2 
->  Ci4Hj^o  +  C2H2 

Irradiation  of  these  product  ions  with  514.5  rm  (2  eV)  photons  results  in  the 


following  photodissociation  processes: 


^6^12*  CgHg'*’  +  CgHg 

->  CieHio^  +  H2 

*-16^10^  *-16^8*  ^2 

^16^11^  ^16^10^  ^ 

The  dimer  ion,  ^16^12^'  photogragment,  C2^gH2^Q''',  react  with 

pheny  lacety  lene  with  rate  constants  of  3.2  x  10”^^  and  3.1  x  10”^^ 
amVmolecule-s,  respectively,  to  form  "trimeric”  ions,  C24Hj^g'''  and  C24H]L6^* 

The  fragment  ions  formed  by  dissociation  of  excited  parent  ions  also 
react  with  pheny lacety lene.  The  following  rate  constants  were  determined: 

3 

on-^/mo  lecu  la-sec 
xlO^O 

Cr^H4‘^  +  CgHgCaCH  ->  C2_2Hg‘^  +  C2H2  12.4  +1.4 

CgH3'^  +  CgHgCaCH  ->  Products  12.2  +1.2 

C6H2'^  +  CgHgCaCl  ->  Products 


21.1  +2.3 


C4H4'''  +  CgHgCsCH  ->  Products  11.4  +1.8 

^4^3^  +  CgHjCsCH  ->  ProdtK:ts  13.9  +0.8 

C4H2'''  +  CgHgCaCH  ->  Products  22.9  +4.3 

That  is,  many  of  the  types  of  reaction  observed  in  our  earlier  study  of  ionic 

reactions  in  diaoetylene  are  observed,  nanely  that  the  highly  unsaturated  ions 
(i.e.  CgH2'*‘  and  C4H2''')  react  at  every  collision  with  the  precursor  acetylenic 

molecule,  and  that  efficient  condensation  reactions  occur  followed  by 
elimination  of  a  molecule  of  acetylene. 

Because  the  details  of  such  ion-molecule  polymerization  and 
polymerization/dissociation  mechanisms  are  largely  controlled  by  thermodynamic 
parameters,  experiments  are  also  carried  out  to  elucidate  the  thermodynamics 
of  the  ionic  mechanisms  in  these  systems.  Proton  transfer  reactions  are 
ubiquitous  in  such  systems,  and  therefore  determinations  of  the  proton 
affinities  of  the  parent  molecules  are  useful  in  understanding  the  overall 
reaction  mechanisms.  The  proton  affinity  of  phenylacetylene  was  determined 
through  "bracketing"  experiments  in  v^ich  a  series  of  bases,  for  vAiich  the 
thermodynamic  parameters  associated  with  proton  donation  are  well  established: 

BH"^  +  CgHgC=CH  ->  +  B 


The  resalts  were  as  follows: 


CgHgCsCH 


i-CgH^yCOOCHg 
(C2Hgj  2^0 
(C2H5) 20 


C2HgCOOCH3 


Proton  affinity (B) 
kcal/mol 

201.6 

201.4 

200.2 

200.2 


Proton  Transfer  Observed? 

Mo 

No 

Yes 

Yes 


The  bracketing  experiments  lead  to  a  value  for  the  proton  affinity  of 
phenylacetylene  of  200.7+0.5  kcal/mol. 
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Phenylethylene  (Styrene) 

Previous  studies  of  ionic  processes  in  styrene  reported  the  occurrence  of 
the  processes: 

CgHs'*’  +  C3HQ  -> 

->  Cj^o*^10*  *  ^6^6 

but  no  kinetic  data  were  reported.  In  the  leinvestigation  of  the  system,  we 
determined  the  rate  constants  and  reaction  mec^ianisnis  for  the  reactions  of 
^8^8^  and  the  fragment  ions  in  order  to  infer  information  about  the  structures 
and  thermochemistry  of  these  fragment  ions  and  their  condensation  products. 
The  following  reactions  were  observed: 

3  "xn 

on  /molecule-sec 
xlO^° 

CgHg  +  CgHgCH*CH2 

-> 

-> 

-> 

CsH?"-  +  CgHgCH-CHj  -> 

(=6^6*  *  CgHgCH^j  -> 

CgHg"^  +  CgH5CH=CH2  -> 

^5^3"  *  C6H5CH-CH2  -> 

^^4^3^  +  CgH5CH=CH2  -> 

Irradiation  of  the  system  with  the  514.8  rin  laser  line  led  to  the  dissociation 
of  the  two  major  product  ions: 

^16^16^  ^8*^8^  *  ^8^8 

^10^10^  ^10^8^  ^2 

Experiments  designed  to  elucidate  the  structures  of  the  and  C2^QHg''' 

ions  were  carried  out.  It  was  reported  in  a  1971  study  of  the  dimerization 
process  in  styrene  that  the  Cj^gHj^g"*"  ion  had  the  structure: 


Products 

‘  Cs"6 

+  C7H7 
Products 
Products 
Products 
Products 
Products 


0.65+0.02 

1.7  +0.2 
15.0  +0.8 

9.1  +1.0 

3.8  +1.1 
13.9  +1.2 
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Thus,  dissociation  of  this  ion  by  loss  of  a  benzene  molecule  might  be  expected 
to  lead  to  the  formation  of  a  species  with  a  structure  such  as; 


In  order  to  test  this  possibility,  the  kinetics  and  photodissociation 
mechanisms  of  the  parent  ions  formed  in  1,2-  and  1,4-dihydronaphthalene  were 
examined;  it  was  concluded  that  these  ions  do  bdiave  like  the  Cj^gHj^g'*’  ion 
generated  in  styrene,  thus  providing  evidence  that  the  structures  given  above 
do  correspond  to  the  species  formed  in  pher^lacetylene. 

The  proton  affinity  of  styrene  was  also  determined  through  "bracketing" 
experiments  like  those  described  above  for  phenylacetylene. 

+  CgHgCH<M2  ->  CgHg"^  +  B 


& 


The  results,  summarized  below,  indicate  that  the  proton  affinity  of  styrene  is 
202.1+0.5  kcal/mol. 


CgH5CH=CH2 


B 

Proton  Affinity 
kcal/mol 

Proton  transfer 
observed? 

(n— 2^ 

202.3 

No 

i-CgHyCOOCHg 

201.6 

Yes 

(C2H5) 200 

201.4 

Yes 

(C2H5)20 

200.2 

Yes 

CoHcCJOOCH 


200.2 


Yes 


structure  and  Reactivity  of  CgHj*  Ions 

Althou^  collision-induced  dissociation  (CID)  studies  have  examined  the 
structures  of  CgHg'*'  ions  formed  as  fragments  in  the  electron  impact  ionization 
of  organic  molecules  and  as  products  of  ion-molecule  reactions,  there  is 
still  conflicting  evidence  about  the  identities  of  these  ions,  whidi  are 
generated  in  the  dissociation  of  the  par^t  ions  of  aromatic  molecules: 


Also,  to  date  no  comprehensive  examination  of  the  kinetics  of  the  reactions  of 
these  ions  has  been  made. 

A  study  was  made  of  the  structures  and  reactivities  of  CgHg'*’  ions  using 
the  ion  cyclotron  resonance  spectrometer.  In  many  systems  (fragment  ions 
formed  in  the  halobenzenes,  or  the  CgHg"^  formed  as  a  result  of  consecutive 
ionic  condensation  reactions  in  acetylene,  for  example)  at  least  two 
populations  of  CgHg"^  can  be  discerned,  one  which  is  highly  reactive  and  one 
whidi  is  relatively  unreactive  with  the  precursor  molecules.  The  identities 
of  these  populations  were  sorted  out  by  generating  CgHg'*’  ions  under  energy 
conditions  such  that  only  the  phenyl  ion  could  be  present,  i.e. 
photodissociation  of  CgHgCl'*'  via  a  two-photon  process  at  514.5  m  must  lead 
exclusively  to  the  formation  of  ions  having  the  phenyl  ion  structure, 
since  there  is  insufficient  aiergy  for  ring  opening  to  occur.  In  this  way  it 
was  demonstrated  that  the  more  reactive  isomer  is  the  phenyl  ion,  contrary  to 
assumptions  made  in  previous  studies.  It  can  be  concluded  that  the  CgHj'*'  ions 
formed  as  products  in  chanical  reactions: 


exhibit  more  acyclic  character  the  higher  the  energy  content  of  the  CgHg'*’  ion 
produced  in  the  reaction. 

The  phenyl  ions  react  with  saturated  hydrocarbons  to  form  a  complex  vAiich 
may  dissociate  to  form  benzene  and/or  protonated  benzene.  Even  throuc^  these 
processes  are  generally  highly  exothermic,  the  reaction  rate  constants  are 
approximately  a  factor  of  ten  lower  than  the  corresponding  ion-molecule 
collision  raite  constants.  Factors  controlling  the  rate  consteuits  and  reaction 
channels  of  ions  interacting  with  organic  and  inorganic  molecules  are 
being  explored  further.  In  experiments  utilizing  partial  deutericn  labelling, 
it  has  been  seen  that  the  phenyl  ion  exhibits  no  selectivity  in  the  hydride 


transfer  reaction: 


CgHg^  +  (CH3)3CD  ->  CgH5{H,D)  +  C4(H,D)g+ 
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MOLECULAR  BASIS  FOR  LASER-INDUCED  VAPORIZATION  OP  REFRACTORY  MATERIALS 


J.W.  Hastle,  D.W.  Bormell,  and  P.K.  Schenck 
High  Temperature  Chemistry  Group 
Ceramics  Division 

Institute  for  Materials  Science  and  Engineering 
National  Bureau  of  Standards 

I.  INTRODUCTION 

The  objective  of  this  project  is  to  develop  a  molecular-level 
understanding  of  the  vaporization  and  related  ionization  and  aggregation 
processes  which  occur  when  high-power  laser  energy  is  dissipated  in 
carbon-containing  materials.  These  highly  complex  processes  can  be  expected 
to  Include  the  thermodynamically  controlled  production  of  C^  (n  s  1  -  10) 
molecular  species,  and  the  non-equilibrium  production  of  higher  cluster 
species  (n  -  10  -  1000),  in  adltion  to  the  formation  of  electrons,  and 
positive  and  negative  ions.  Inorganic  carbon-based  materials,  such  as 
graphite,  metal  carbides,  and  carbon-carbon  composites,  have  intrinsically 
high  resistance  to  thermal  degradation.  This  behavior  results  from  an 
inherently  high  sublimation  energy,  high  melting  point,  and  high  thermal 
conductivity.  The  development  of  effective  laser-resistant  materials  and 
obscurants  is  expected  to  rely  mainly  upon  such  carbon-based  materials  [Baker 
and  Hager,  1984;  Stlckley,  1985]. 

Despite  considerable  research  over  the  past  several  decades,  the  basic 
thermodynamic,  kinetic,  and  mechanistic  description  of  the  thermal 
decomposition  of  graphite  and  related  materials  is  still  in  considerable 
doubt.  In  fact,  two  widely  used  but  discrepant  sets  of  thermodynamic  tables 
exist  for  the  carbon  vapor  species,  the  JANAF  [1971-82]  thermochemical 
tables  and  the  so-called  Livermore  tables  [cf.  Baker,  et  al,  1983].  Even  the 
melting  point  of  graphite  is  a  subject  of  considerable  controversy. 
Experiments  have  been  conducted  on  the  effect  of  moderate  power  CO2  and 


Nd-YAG  lasers  on  the  vaporization  of  various  foras  of  graphite  [Covington,  et 
al,  1977;  Kirillin,  et  al,  1984;  Lincoln  and  Covington,  1975;  Lundell  and 
Dickey,  1977;  Lundell,  1982;  Meyer,  et  al,1973;  Milne,  et  al,  1972; 
Schelndlin,  1984;  Zavitsanos,  1968].  Most  of  these  studies  have  concentrated 
on  the  macroscopic,  non**aiolecular  processes  such  as  crater  formation, 
material  removal  rate,  plume  dynamics,  and  phase  equilibria.  We.  have 
reviewed  the  recent  literature  [Hastle,  et  al,  1984]  and  find  that  there  Is 
considerable  disparity,  particularly  concerning  relative  and  absolute  partial 
pressures  of  the  various  carbon  polymer  species. 

The  work  reported  here  has  utilized  a  very  high  pumping  speed  mass 
spectrometer  system,  with  a  Nd/YAG  laser  providing  the  vaporizing  heat 
source.  Special  emphasis  has  been  given  to  the  relationship  of  the  C^ 
species  partial  pressures  to  laser  parameters  such  as  fluence,  Irradiation 
angle,  and  Impact  repetition. 

II.  EXPERIMENTAL  APPARATUS  AND  PROCEDURE 

Installation  of  a  new  20  Hz  Nd~YAG  laser  system,  and  its  coupling  to  the 
high  pressure  sampling  Mass  Spectrometer  (MS),  was  recently  completed.  The 
experimental  setup  is  otherwise  similar  to  that  reported  previously  [Hastie, 
et  al,  1984],  and  is  shown  schematically  in  figure  1.  The  salient  feature  of 
the  sample  mount  region  is  the  very  high  pumping  speed,  exceeding  2000  l/sec 
at  pressure  up  to  10“5  atm  (1  atm  =  101325  Pa). 

The  laser  system  operates  at  wavelengths  of  355,  532,  or  1064  nm  with 
energies  up  to  200  mJ  per  pulse  at  1064  nm.  The  present  studies  utilized 
532  nm  as  this  wavelength  provided  the  best  match  of  energy  and  vaporization 
conditions  (10-40  mJ  into  a  250  urn  spot).  For  physical  access,  the  beam  is 
perlscoped  to  the  level  of  the  sample  and  then  piped  to  a  motor-driven  500  mm 
focusing  lens  via  dielectric  mirrors  which  are  housed  in  a  nitrogen-purged 


9  X  12  cm  conduit  (not  shown  In  figure  1).  A  radiometer  was  used  to  directly 
measure  the  laser  fluence,  characterize  the  transmission  losses  of  the  beam 
through  the  conduit  system,  and  calibrate  lamp  power  to  delivered  laser 
power.  Radiometric  measurements  were  made  periodically  during  experimental 
runs  to  confirm  constant  laser  power  conditions. 

Alignment  of  the  laser  impact  point  with  the  mass  spectrometer 
beam-acceptance  axis  was  monitored  using  a  photodiode  located  Just  above  the 
mass  filter  ion  source.  Optimum  alignment  was  obtained  by  adjusting  the 
final  laser  mirror  to  maximize  the  Intensity  of  the  reflected  laser  light 
from  the  sample  surface.  This  technique  was  shown,  using  a  He/Ne  laser  and 
mirror,  to  provide  accurate  positioning  of  the  laser  impact  point  at  the 
center  of  the  molecular  beam  acceptance  cone  of  the  MS  system.  Because  of 
the  extreme  forward-peaking  of  the  emitted  molecular  beam,  the  acceptance 
angle  of  the  beam  system  was  set  at  about  13*,  with  a  usable  targeting  point 
radius  about  the  nominal  beam  center  origin  of  -0.3  cm.  From  an  inspection 
of  the  deposit  patterns  on  the  beam  wall  directly  above  the  sample  point,  the 
beam  distribution  appeared  to  be  much  more  collimated  than  COS^  (typical  of 
effusion)  due  to  the  beam  self-focusing.  The  beam  centerline  emission  angle 
was  observed  to  be  within  -±  2*  of  normal  to  the  surface  with  the  sample 
angle  at  15“  from  the  horizontal  laser  beam.  Both  forward-peaking  and  normal 
beam  emission  were  independent  of  laser-sample  angle  and  of  whether 
individual  laser  shots  always  struck  a  fresh  surface,  or  up  to  18,000  shots 
were  allowed  to  impact  the  same  spot. 

Most  of  the  recent  experiments  used  high  density/high  purity  graphite  of 
similar  overall  purity  to  the  spectroscopic  grade  graphite  used  previously. 
The  primary  reason  for  changing  the  specimen  source  was  to  obtain  a  larger 
sample  cross  section  (-5.6  cm).  Specimens  0.24  cm  thick  were  parted  from  a 


large  bar  with  a  carefully  cleaned  steel  cutting  tool.  The  polishing  step 
was  omitted  to  eliminate  possible  contamination  by  silicon  carbide,  even 
though  this  did  not  appear  to  be  a  problem  in  our  earlier  studies.  The 
surface  roughness  of  the  present  samples  was  only  slightly  greater  than  that 
of  the  ll^tly  polished  spectroscopic  grade  material  used  earlier.  A  mass 
spectral  survey  Indicated  no  contaminant  peaks  other  than  those  noted  for  the 
spectroscopic  rod  source  (l.e.,  only  minor  hydrocarbon  content). 

The  experimental  data  acquisition  process  consisted  of  the  following 
steps: 

1.  Position  sample  to  the  limit  position  of  the  translator,  adjust  beam 
position  and  focus. 

2.  Position  sample  to  location  of  Interest. 

3.  Tune  mass  filter  to  mass  of  Interest. 

4.  Open  laser  shutter,  letting  pulses  strike  surface. 

5.  Collect  ion  signal  during  the  -  5ms  time  period  which  follows  each 
pulse  using  a  signal  averager  with  1024  data  sampling  points. 

Between  100  and  1000  shots  were  generally  averaged  to  obtain  good 
3 1 gnal-to-nol se . 

6.  (a)  For  the  static  mode,  successive  laser  shots  Impact  the  same 

location. 

(b)  For  the  dynamic  mode,  the  sample  Is  moved  after  every  shot. 

The  maximum  movement  was  about  lOum  between  shots. 

7.  Repeat  steps  2-6  for  each  mass  position  of  interest. 

This  procedure  produces  mass  selected  profiles  of  ion  intensity  versus  time, 
as  shown  by  the  example  given  in  figure  1.  The  time  scale  includes  the 
laser  on/cff  pulse  (10  ns),  together  with  the  much  longer  (millisecond)  times 
for  sample  cooling  and  beam  transit  to  the  mass  spectrometer.  Earlier  work 


CHaatie  et  al,  1984]  indicated  that  these  profiles  contain  data  on  neutral 
species  flight  times  and  hence  post-expansion  vapor  temperature,  in  addition 
to  the  surface  thermal  history.  However,  the  earlier  use  of  multiple  laser 
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shots  at  each  location  produced  species  Intensities  which  apparently 
corresponded  to  unit  accomodation  coefficient  due  to  the  formation  of  a 
crater  which  could  behave  as  a  crude  Knudsen  cell.  With  successive  laser 
shots  at  a  fixed  location,  signal  Intensities  decreased  roughly  in  proportion 
to  the  logarithm  of  the  number  of  shots.  We  interpreted  this  effect  as  being 
due  to  crater  formation  causing  an  enhanced  forward-beam  peaking,  with  a 
consequent  loss  of  signal  off  the  beam  centerline. 

During  this  reporting  period,  we  have  devised  and  installed  a  two  axis 
translation  system.  Operating  under  real-time  control  with  in-house 
developed  code,  this  system  allows  samples  to  be  positioned  or  rastered,  as 
desired,  between  laser  shots.  The  system  permits  positioning  over  a  square 
area  2.5  cm  on  a  side.  Raster  rates  up  to  200  um/sec  are  possible, 
corresponding  to  somewhat  over  90  percent  shot-to-shot  overlap  with  a  -250ii 
beam  spot  dimension.  With  this  new  raster  system,  it  is  now  possible  to 
monitor  the  vaporization  process  continuously,  from  the  pristine  surface 
condition  to  the  cratered  condition  of  our  earlier  work. 

III.  RESULTS  AND  DISCUSSION 

Experiments  were  carried  out  during  this  reporting  period  to  examine  in 
greater  detail  the  special  features  noted  earlier  for  the  vaporization  of 
graphite.  In  particular,  the  observation  of  thermodynamic  equilibrium  at 
nearly  unit  accomodation  coefficient  from  cratered  surfaces  was  tested  by 
operating  in  the  raster  mode,  where  each  shot  was  directed  at  a  fresh 
surface.  The  earlier  observation  of  anomalously  high  C^/Cj  ratios, 
attributed  to  non-thermal  sources,  was  also  reexamined. 


III.1  Scanning  Veraus  Static  Laaer-Surfaee  Interaction 

Our  Initial  mass  spectral  observations  showed  a  decreasing  vaporization 
rate  with  increasing  number  of  laser  shots  at  a  fixed  surface  location.  The 
effect  of  laser-induced  vaporization  from  a  fresh,  non-cratered  surface  was 
examined  by  scanning  the  laser  beam  across  the  sample  surface.  Thus  each 
sample  location  received  only  a  small  (-20-25)  number  of  laser  shots. 

Figure  2  shows  the  effect  of  systematically  varying  the  number  of  laser  shots 
on  the  mass  spectral  Ion  Intensities  for  C*  (n  -  1  -5).  Note  that  the  rate 
of  decrease  in  signal  (normalized  with  respect  to  the  number  of  laser  shots) 
with  number  of  laser  pulses  Is  considerably  more  pronounced  for  C|  and  Cg 
than  for  and  C^.  The  higher  mass  species  intensltes  are  essentially 
constant,  within  error  limits,  particularly  after  only  a  few  thousand  shots. 
These  observations  suggest  that  during  the  early  stages  of  vaporization, 
under  continuous  laser  Irradiation,  the  lower  mass  species  may  indeed 
predominate.  This  effect  could  account  for  the  disparate  literature  results 
on  the  relative  abundance  of  the  various  C^j  species  (e.g.,  see  Clarke  and 
Fox,  [1969];  Wachl  and  Gilmartln  [1970]).  Wachi  and  Gilmartin  [1970] 
reported  a  similar  effect,  with  Increasing  strongly  relative  to  and  C2 
over  time,  which  they  attributed  to  changes  In  surface  morphology  and/or  the 
possible  formation  of  a  carbon  allotrope.  Section  III. 5  discusses  our 
surface  analysis  results. 

Table  1  shows  a  comparison  of  the  present  laser-scanned  results  with  the 
earlier  [Hastie,  et  al,  1984]  static  data.  Note  that  the  species  ratios  are 
in  reasonable  agreement  between  the  two  cases,  except  for  and  the  reason 
for  this  la  currently  under  Investigation  (see  section  III. 5).  The  slightly 
lower  value  for  is  unexpected  and  may  result  from  an  Instrumental  effect 
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on  peak  widths.  This  peak  width  reproducibility  problem  is  under  investiga¬ 


tion,  but  is  probably  caused  by  insufficient  bandwidth  in  the  data  collection 
process  (see  section  III. 3).  The  direct  observation  of  Cg  -  Cg  reflects  the 
improvement  in  sensitivity  we  have  accomplished  with  the  new  experimental 
arrangement.  Cg  currently  represents  our  detection  limits,  with  slgnal/nolse 
ratios  of  about  1 .  A  significant  Improvement  in  signal-to-nolse  can  be 
expected  when  we  convert  to  signal  counting  detection  techniques  (see  section 
III. 3  below). 

111. 2  Evidence  for  Local  Thermodynamic  Equilibrium 

In  the  previous  static  experiments  [Hastie,  et  al,  1984],  it  was  noted 
that  the  various  carbon  polymer  species  appeared  to  be  generated  with  a 
common  source  (sample  surface)  temperature,  T.  Regardless  of  the  exact  model 
for  the  vaporization  and ■ beam- forming  process.  If  a  collection  of  gas 
molecules  are  in  thermal  equilibrium,  their  velocity  function  is  proportional 
to  (Tq/m)''^^,  where  M  is  the  individual  species  molecular  weight  and  is 
the  gas  temperature.  Hence  a  plot  of  time-of-arrival,  At  (which  is  inversely 
proportional  to  velocity)  versus  (M)'^^  should  show  a  linear  correlation  for 
thermal  equilibrium  conditions.  Figure  3  shows  the  peak  (most  probable) 
arrival  time  of  carbon  polymer  species  up  to  for  three  ranges  of 
accumulated  laser  shot  conditions.  Within  repeatability,  the  expected  linear 
relation  is  found.  The  outlier  of  Cj,  at  14,000  shots  in  this  particular  data 
set  also  showed  an  anomalously  wide  intensity-time  profile  which  was  not 
reproduced  in  other  experimental  series. 

111. 3  Relative  Temperatures  and  Instrumental  Response 

In  order  to  obtain  good  time  resolution,  the  signal  averaging  has  been 
performed  by  collecting  the  amplified  multiplier  output  current  (analog)  of 
the  MS.  The  analog  signal  is  then  digitized  by  a  fast  8  bit  digitizer  in  the 


fflultichannel  analyzer  (MCA)  and  accumulated  In  individual  channels.  This 
process  results  in  imposing  the  final  DC  amplifier  risetime  function  (with  a 
time-constant  of  the  order  of  a  few  kHz)  on  the  intensity  profiles,  as  well 
as  subjecting  low  level  signals  to  "noise”  caused  by  digitizer  granularity. 

A  more  desirable  method  would  be  to  use  the  electron  multiplier  in  the 
saturated  mode,  and  count  individual  ion  events.  The  currently  available  MCA 
has  unacceptable  dead  time  requirements  in  its  multichannel  scaler  mode  at 
minimum  channel  dwell  time  (-5  ya/channel).  It  should  be  noted  that  the 
channel  time  is  already  a  compromise  and  should  be  faster,  by  about  an  order  of 
magnitude,  for  optimum  resolution  of  the  temperature  information. 

To  Improve  the  signal  detection  time  constant,  we  decreased  the  DC  gain 
by  a  factor  of  100,  by  changing  the  amplifier  load  resistor  from  10^  0  to 
10^  Q,  which  Increased  the  system  bandwidth  to  -10^  Hz.  Preliminary  results 
indicated  that  with  normal  amplification  of  the  DC  signal,  risetime 
electronic  delay  effects  are  very  significant.  The  reduction  in  integrated 
signal  intensity  was  nearly  equal  to  the  factor  of  100  reduction  in  system 
gain.  However,  peak  widths  were  broadened,  possibly  causing  excessive 
Integrated  signal  reproducibility  problems.  For  example,  the  width  of  the 
12  amu  signal,  inherently  the  most  narrow  ion  intensity-time  profile,  was 
almost  a  factor  of  two  narrower  at  the  lower  gain  setting.  In  addition,  the 
most  probable  arrival  time  for  C.j  was  100  us  less.  Both  these  results 
indicate  at  least  a  100  us  characteristic  (instrumental)  delay  for  the  higher 
gain  setting.  Operation  at  the  lower  gain  reduces  system  sensitivity  such 
that  the  C5  species  data  are  nearly  lost  in  the  noise.  Thus,  before  absolute 
time-of-arrival  data  can  be  collected,  and  detailed  modeling  of  the  time- 
resolved  intensity  profile  data  can  be  performed,  pulse-counting  using 
multichannel  scaling  will  have  to  be  implemented  to  preserve  system 


sensitivity  and  eliminate  the  system  response  delay.  Consequently,  although 
relative  temperatures  derived  by  various  models  from  the  low  bandwidth  data 


are  valid,  absolute  temperatures  are  not  easily  extracted  from  the  present 
results. 

III. 4  Effect  of  Laser  Irradiation  Energy  (Fluence) 

A  series  of  experiments  were  conducted  to  quantitatively  determine  the 
effect  of  laser  fluence  on  the  mass  spectrum  of  carbon  polymer  species.  It 
was  noted  earlier  [Hastle,  et  al,  1984]  that  different  experimental  series, 
with  only  approximately  equal  laser  irradiation  energy,  all  yielded  the  same 
thermodynamic  temperature,  4200  ±  300  K.  This  effect  could  be  due  to  a 
physicochemical  constraint,  such  as  a  solid-liquid  phase  change.  Graphite 
melting  temperatures  are  believed,  from  literature  studies,  to  fall  in  the 
range  4000-4500  K, 

The  effect  of  a  three-fold  variation  of  laser  fluence  on  the  arrival 
time  to  molecular  weight  dependence  is  given  in  figure  4.  From  the  curve 
slopes  in  figure  4,  the  temperature  is  calculated  to  be  essentially  constant 
(see  section  IV).  Since  peak  widths  also  did  not  vary  with  laser  power,  the 
total  cooling  time  also  did  not  apparently  change  signf icantly.  A  simple  T^ 
energy  loss  argument  shows  that  if  only  the  same  area  were  heated,  the 
temperature  should  have  varied  by  a  factor  of  1.3,  e.g.  from  4200  to  5500  K, 
if  4200  K  represents  the  minimum  temperature.  Mote  that  the  laser  beam  focus 
was  carefully  maintained  throughout  the  series,  and  the  surface  was  scanned 
to  eliminate  possible  cratering.  Thus  we  infer  that  the  heated  area  changes 
directly  with  laser  fluence  at  a  temperature  fixed  by  the  system  constraint. 
Considering  the  short  time  scale,  the  thermal  transport  probably  occurs  in  a 


very  narrow  surface  zone.  Otherwise  we  would  expect  to  see  a  dependence  of 
tliBe-of-arrlval  on  laser  fluence.  The  observation  that  the  C2  -  C5  species 
intensities  vary  directly  with  fluence  also  supports  this  argument. 

The  intensity  showed  an  unexpected  square  law  dependence  on  the  laser 
fluence,  an  indication  that  is  also  produced  by  a  photodissoclation  or 
similar  non-thermal  mechanism.  This  result  is  particularly  surprising  in 
view  of  our  special  efforts  to  minimize  laser  interaction  with  the  plume — the 
laser  Impacts  at  a  shallow  angle  (15^),  the  laser  pulse  time  is  very  short 
(-10  nsec),  and  the  laser  Is  relatively  tightly  focused  (-250  u  diameter  at 
Impact).  Hence  the  photodissoclation  process  must  be  occuring  very  early  In 
the  vaporization  process,  or  be  a  laser/surface  Interaction  phenomenon.  From 
simple  gas  velocity  arguments,  the  known  maximum  surface  temperature  of 
4200  K  [Hastle,  et  al  1984],  and  the  10  nsec  laser  on-time,  the  entire 
laser/gas  interaction  must  occur  at  a  distance  within  30  pm  of  the  surface. 
The  observation  that  peak  width  is  also  insensitive  to  laser  fluence  implies 
either  that  the  laser  heated  spot  cools  at  a  rate  nearly  as  fast  as  it  heats, 
or  that  the  photoproduction  mechanism  is  a  result  of  the  laser/surface 
interaction  and  persists  throughout  the  cooling  time.  This  process 
undoubtably  is  responsible  for  some  of  the  earlier  mentioned  dlscrepencles  in 
the  literature  concerning  species  ratios.  Also,  the  practical  requirement  of 
an  enhanced  monomer  production,  to  maximize  energy  dissipation  in 
laser-resistant  shielding,  may  be  met  naturally  as  a  consequence  of  the 
extremely  high  fluence  levels  likely  to  be  present. 


III. 5  Surface  Characterization 


The  apparent  change  in  the  vaporization  kinetics  and  beam  formation 
process,  as  a  function  of  the  number  of  laser  pulses,  prompted  a  study  of 
the  surface  morphology.  Both  scanning  electron  microscopy  (SEM)  and  surface 
profilometry  were  employed  to  study  the  laser  irradiated  samples. 

In  typical  laser-scanning-mode  experimental  runs,  the  samples  were 
translated  at  -200  pm/sec,  which  yields  about  95  percent  overlap  of  each 
successive  shot.  A  portion  of  the  sample  was  irradiated  under  these 
conditions  in  a  single  pass.  Preliminary  SEM  observation  of  the  surface 
indicated  that  the  lack  of  fiducial  marks  made  it  difficult  to  even  Identify 
the  path  of  the  laser.  Figure  5  is  an  SEM  photograph  of  a  portion  of  the 
laser  path  identified  by  two  fiducial  marks  manually  scribed  approximately 
perpendicular  to  the  laser  path,  which  was  readily  visible  to  the  naked  eye 
by  an  interference  produced  color  change  on  the  sample  surface.  The  15* 
angle  of  Irradiation  should  produce  an  -400-500  pm  "footprint"  on  the  sample. 
With  the  aid  of  the  manually  scribed  fiducial  marks,  such  a  "footprint"  is 
barely  observable  in  the  SEM  micrograph.  The  graphite  grain  size  can  be  seen 
to  be  on  the  order  of  100  pm  in  the  sample.  This  small  surface  change  is 
quite  different  from  the  craters  generated  in  the  earlier  static  experiments. 
Such  craters  were  typically  sharp-edged  pits  of  considerable  depth  relative 
to  the  spot  size. 

The  combination  of  surface  roughness  and  low  contrast  in  the  SEM 
micrograph  of  figure  5  makes  it  impossible  to  estimate  the  rate  of  sample 
removal  for  the  mass  spectrometer  experiments.  In  an  attempt  to  provide  such 
data,  an  automated  profllometer  study  was  made  of  the  surface  morphology 
of  a  portion  of  the  single  pass  laser  path  on  the  sample.  The  profllometer 
made  repeated  scans  of  the  sample  across  the  laser  path  at  10  pm  intervals 
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for  a  distance  of  4  on.  The  result  of  100  scans  Is  shown  In  figure  6.  This 
3*0  Isometric  plot  of  the  resultant  data  covers  a  1  by  4  mm  portion  of  the 
surface.  The  vertical  scale  la  about  75  iim  from  the  bottom  of  the  slice  to 
the  average  of  the  surface  peaks.  From  this  plot,  the  surface  roughness  can 
be  seen  to  be  on  the  scale  of  100  urn  across  the  surface  and  25~50  pm  into  the 
surface.  The  laser  path  is  not  discernable. 

IV.  MODELING  OF  THE  THERMAL  VAPORIZATION  PROCESS. 

The  original  velocity  analysis  used  to  establish  the  equilibrium  nature 
of  the  vaporization  process  of  the  vaporized  neutral  species  is  an  extremely 
simplified  model.  As  we  noted  earlier  [Hastle,  et  al,  1984],  the  slope  of  At 
versus  provides  a  measure  of  temperature,  but  is  unreasonably  low. 

Previous  analysis  of  the  data  equated  a  maximum  beam  temperature  with 
the  fastest  arriving  species.  A  more  reasonable  method  to  model  the 
gasdynamlcally  expanded  beam  temperature  is  to  deduce  the  temperature  from 
the  time-of-f light  data  by  assuming  that  the  beam  species  obey  a  Maxwellian 
velocity  distribution.  The  characteristic  velocity  of  a  Maxwellian  velocity 
distribution  is  the  most  probable  velocity  as  given  by: 

v„p  -  (2kT/M)’/2 

where  k  is  Boltzraan  constant,  T  is  temperature,  and  M  is  mass.  Converting  the 
velocity  distribution  to  an  arrival  time  distribution,  by  assuming  Instantaneous 
source  and  response  functions,  gives  a  most  probable  arrival  time  of  : 

t^p  -  L 

where  L  is  the  path  length  from  the  sample  to  the  mass  analyser.  As  in  the 
previous  analysis,  the  assumption  of  thermal  equilibrium  requires  that 
t^p  vs.  /m  is  linear.  The  beam  temperature  can  be  deduced  from  the  slope  of 
a  line  drawn  through  the  points  of  graphs  in  figures  3  and  4  and  la 
proportional  to  the  slope  squared.  Only  the  high  bandwidth  data  (10^  Q) 


reasonably  approximate  the  instantaneous  response  assumption.  Analysis  of 
the  data  of  figure  4  results  in  a  temperature,  averaged  over  all  the  data,  of 
290  ±  50  K.  This  beam  temperature  is  equivalent  to  a  supersonic  expansion 
cooling  factor  of  about  14  (based  on  our  surface  temperature  of  4200  K).  This 
degree  of  cooling  indicates  an  essentially  complete  gasdynamic  cooling  from 
an  extremely  small  region.  Such  a  conclusion  is  in  accord  with  the  small 
focus  zone  ('250  pm)  and  high  total  graphite  vapor  pressure  ('1  atm). 

The  time-resolved  intensity  profiles  actually  contain  the  entire  thermal 
history  of  the  laser  vaporization  process.  A  better  basis  for  recovering 
that  information  la  possible  by  deconvolution  of  the  entire  time-of-flight 
signal,  which  la  an  Integral  of  a  vapor  generation  function,  a  velocity 
distribution  function,  and  the  mass  analyser  detector  function; 

t  6 
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The  functions  in  this  convolution  are; 

g(t)  is  the  Instrument  detector  and  electronics  response  function; 

f(t)  is  the  arrival  time  distribution  function  e.g.  for  a  Maxwellian 
velocity  distribution; 

C(T)  is  the  species  creation  function,  dependent  on  temperature, (T)  and 
time  (t); 

K  Ja  the  calibration  factor  for  the  number  density  detector; 

I(t)  is  the  observed  time  dependent  ion  current  signal. 


Such  a  deconvolution  has  been  treated  by  U.S.  Young  [1975]  for  the 
simple  case  of  a  mechanically  chopped  isothermal  beam.  He  showed  that  the 
above  convolution  accurately  models  the  obs^ved  mass  spectrometer  signals. 
Current  efforts  are  directed  to  the  development  of  a  parametric  comput«'^  code 
to  solve  the  deconvolution  by  numerical  integration  for  the  more  complex  case 
of  a  tlme~varying  (exponential  decay)  beam  source  temperature,  and  a 
Clauslus-Clapeyron  type  species  generation  function  in  place  of  the  simple 
on-off  modulation  function  of  Young  [19753.  As  a  parametric  model,  it  can 
allow  for  eventual  inclusion  of  explicit  functions  for  macroscopically 
Important  effects  such  as  substrate  thermal  conductivity,  radiation  balance, 
and  thermal  transport  by  material  loss.  However,  the  surprisingly  large 
influence  of  the  analog  MS  detector  response  function  (see  III. 3)  will  have 
to  be  drastically  reduced  before  the  model  can  be  implemented.  The  strong 
effects  of  the  current  response  function  act  as  a  large  noise  background  to 
the  convolution,  making  the  model's  sensitivity  to  important  experimental 
variables,  such  as  cooling  rate,  too  poor  to  yield  reliable  results  at 
present.  Methods  to  minimize  the  instrument  response  function  using  counting 
techniques  are  being  investigated. 

V.  SUMMARY  AND  PLANS 


'/.  1  Summary 

A  new  laser  system  has  been  installed  and  tested.  In  addition,  a 
two-axis  sample  stage,  which  allows  for  single  and  multiple  shot  measurements 
of  selected  areas  of  the  sample's  surface,  has  been  installed  and 
demonstrated.  The  laser  to  sample  to  mass  spectrometer  geometry  is  similar 
to  that  of  our  previous  work,  where  the  laser  impacts  the  surface  at  a  shallow 
angle  (-15»)  to  minimize  laser/plume  interactions.  Emission  of  vapor 


species,  however,  was  noted  to  occur  perpendicular  to  the  surface,  was  highly 
forward-peaked,  with  a  much  narrower  colllmatlon  than  a  COS^  distribution. 
These  effects  were  noted  even  from  an  uncratered  surface. 

The  observation  of  thermal  equilibrium,  reported  by  us  previously  for 
laser-pitted  surfaces,  was  also  shown  to  be  characteristic  of  a  surface  with 
"natural"  roughness  (-25  ub).  Abnormally  high  concentrations  of  and,  to  a 
lesser  extent,  C2,  relative  to  the  other  species,  were  found  with  fresh 
surfaces.  This  result  is  perhaps  one  cause  of  the  disparate  literature 
results  for  the  relative  partial  pressures  of  carbon  species.  Another 
contributing  factor  In  this  regard  is  the  square-law  dependence  of 
Intensity  on  laser  fluence.  This  behavior  Is  tentatively  attributed  to  a 
photodlssoclatlon  process  occurlng  on  the  surface.  It  is  difficult  to 
distinguish  between  a  laser-vapor  Interaction  and  a  persistent  surface  effect 
until  we  measure  the  characteristic  cooling  time  of  the  surface.  The  role  of 
surface  roughness  la  clearly  a  significant  factor  which  needs  further 
Investigation. 

The  temperature  reached  by  the  graphite  surface,  of  >>200  ±  300  K,  is 
relatively  Independent  of  laser  fluence  over  at  least  a  three-fold  power 
range  at  532  nm.  It  appears  that  the  temperature  could  be  constrained  by  a 
physicochemical  process,  e.g.  the  thermal  halt  due  to  a  solid-liquid  phase 
change.  This  result  also  indicates  that  the  kinetics  of  the  beam  expansion 
and  the  internal  temperature  of  the  plume  are  Independent  of  laser  energy. 
Experiments  testing  the  contribution  of  the  mass  spectrometer  detector 
response  function  to  the  shape  of  the  tlme-of-arrival  spectra  show  that,  at 
instrument  gains  necessary  to  good  sensitivity,  the  bandwidth  limitations  are 


of  the  same  order  as  the  tenperature  effects.  Thus  the  tloe-of-arrival  data 
currently  provide  relative  temperatures,  with  the  prospect  of  extracting 
absolute  temperatures  yet  to  be  resolved. 

V.2  Plans 

In  FY'-SS,  we  intend  to  continue  studies  on  the  laser  energy  dependence 
of  the  time  amd  mass  resolved  molecular  species  from  graphites.  'In  order  to 
decouple  the  instrument  response  function  from  the  tlme-of-ar rival  data,  the 
MS  detector  system  will  be  adapted  to  pulse  counting,  using  an  appropriate 
multichannel  scaler.  This  will  allow  for  model  analysis  of  time~of*arrlval 
data  to  identify  and  extract,  if  possible,  the  source  temperature  from  the 
expansion  Maxwellian  velocity  distribution. 

Preparation  will  begin  for  adapting  the  vacuum  system  for  ion  detection. 
This  will  Involve  designing  the'  ion  optics  to  allow  passage  of  thermally 
ionized  species  through  the  differentially  pumped  apertures. 

The  surface  characterization  results  indicate  that  the  sample 
preparation  procedure  (see  section  II)  Is  clearly  not  adequate  for  these 
studies  if  crater  formation  and  material  removal  rate  information  is  to  be 
obtained.  For  this  reason,  diamond  polished  planchets  of  high-purity 
graphite  have  been  obtained.  The  mirror  finish  on  these  samples  will  be 
characterized  by  both  SEM  and  profilometry  prior  to  laser  irradiation. 
Preliminary  SEM  observations  indicate  no  structural  features  above  submicron 
size.  By  eliminating  the  surface  roughness  effects  Jn  the  vaporization 
process,  it  will  be  possible  to  measure  accurate  mass  removal  rates,  and 
provide  the  gravimetric  calibration  necessary  to  obtain  accurate  partial 
pressures  from  scanned  surfaces.  It  should  also  be  informative  to  follow  the 
cau’bon  species  distributions  as  vaporization  proceeds  from  a  clearly  surface 
dominated  region  to  the  cratered  condition  of  the  static  experiment  samples. 


As  surface  roughness  appears  to  play  a  significant  role  in  species 
distributions,  particularly  the  low  oass  species  of  importance  in  practical 
applications,  graphites  with  submicron  surface  smoothness  will  be  studied.  I 
order  to  investigate  the  role  of  surface  finish,  specimens  will  be  examined 
after  laser  irradiation  to  determine  where  the  vaporization  occurs  (perhaps 
at  grain  boundaries) . 

The  slight  off-axis  sampling  geometry  used  (see  section  II)  may  favor 
beam  segregation  with  an  apparent  loss  of  the  heavier  species.  For  this 
reason,  we  are  modifying  the  MS  system  to  allow  for  true  perpendicular 
(on-axls)  molecular  beam  sampling.  We  Intend  to  test  geometries  where  the 
laser  beam  is  also  perpendicular  to  the  sample  surface,  and  hence  co-axial 
with  the  molecular  beam,  in  addition  to  conditions  where  the  beam  strikes  at 
low  angles  to  a  surface  perpendicular  to  the  MS  beam  axis.  We  will  also 
examine  the  possibility  of  folding  the  laser  path  within  the  vacuum  system 
for  perpendicular  beam  extraction  and  low  angle  laser  irradiation. 

Because  of  the  difficulty  of  real-time  pyrometry  on  a  submicrosecond 
time  scale,  we  intend  to  investigate  the  persistence  time  of  the  residual  hot 
spot  on  the  surface  after  the  laser  pulse  terminates.  Filtered  high  speed 
photodiodes  will  be  used  to  used  to  observe  the  hot  spot  and  the  plume,  with 
simultaneous  recording  of  photodiode  output  and  MS  signal  to  allow 
correlation  of  the  vaporization  process  with  the  laser  pulse  and  subsequent 
history. 

Because  metal  carbides  are  a  attractive  potential  protective  material, 
initial  survey  studies  of  TaC,  as  an  example,  will  also  be  undertaken. 
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d)  Error  estimate  based  on  similarity  to  C_  and  C-  uncertainty. 


Figure  Captions: 


(1)  Schematic  of  laser-lnduced  vaporization  mass  spectrometrlc  facility 
(side  view),  showing  laser-sample^molecular  beam  geometries.  The 
microcomputer  system  shown  controls  the  signal  averager  and  the 
instantaneous  position  of  the  sample  mount/translator.  Completed 
scans  are  recorded  and  then  replotted  by  a  minicomputer.  A  typical 
recent  tlme~resolved  Intensity  profile  for  mass  12  amu  is  shown. 
Indicating  good  slgnal/nolse. 

(2)  Dependence  of  Integrated  signal  intensity  (per  laser  shot)  versus 
accumulated  number  of  laser  shots  for  a  "static"  (see  text) 
graphite  sample  position.  Data  were  recorded  for  Cjj  (n«1  to  5) 
for  1000  shots  each,  then  repeated  to  monitor  the  fall-off  In 
signal  intensity  with  accumulated  number  of  laser  shots.  The  laser 
fluence  was  held  constant  at  37  mJ/shot  during  the  runs. 

(3)  Dependence  of  the  most  probable  tlme-of -arrival  (signal  peak) 
versus  the  square  root  of  mass  for  C^  (n-1  to  5)  ions.  The  linear 
correlation  suggests  local  thermodynamic  equilibrium  of  the  beam 
species.  This  equilibrium  is  essentially  unchanged  with 
accumulated  number  of  laser  shots  on  a  static  graphite  sample.  A 
hlgh-galn  (10^  Q)  input  impedance  was  used. 

(4)  Dependence  of  the  moat  probable  time-of-arrival  (signal  peak) 
versus  the  square  root  of  mass  for  (n-1  to  5)  ions  for  varying 
laser  fluence  on  a  scanned  graphite  sample.  The  input  gain 
impedance  was  10^  q.  Note  that  the  variation  of  arrival  time  with 
laser  fluence  is  less  than  the  repeatability  demonstrated  by  the 
replicate  set  of  37  mJ  data. 

(5)  SEM  micrograph  of  laser  track  across  a  scanned  graphite  target.  The 
laser  path  (see  arrow)  Is  approximately  perpendicular  to  the  two 
deep  manually  scribed  fiducial  marks  which  are  the  most  prominant 
features  in  the  micrograph.  The  two  arrows  mark  the  direction  and 
apparent  edges  of  the  laser  track.  The  scale  marker  is  100  um. 

(6)  Isometric  plot  of  data  obtained  from  an  automated  profilometer 
study  of  the  laser  path  on  the  scanned  graphite  target  of  figure  5. 
The  base  of  the  plot  is  1  by  4  mm  and  the  average  vertical  height 
is  75  um.  The  laser  path  is  parallel  to  the  1  mm  side,  as  shown  by 
the  arrow.  However  the  path  cannot  be  discerned  in  this  figure. 

See  text  for  discussion. 
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